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Introduction 


In a recent article on the mechanical behavior of continuous media [/], the 
physical concept of a fluid was given a mathematical definition. We believe that 
this definition covers almost all real fluids (whenever thermal and other non- 
mechanical effects can be disregarded) and is more general than any proposed 
previously because it accounts for all hereditary effects including stress relaxation. 
The general fluids covered by this definition may even exhibit, in some situations, 
physical phenomena which are usually attributed to solids. 

Here we show that certain steady flow problems can be solved assuming 
only the definition of a general fluid and incompressibility. The class of problems 
considered includes most of the classical flows: simple shearing flow, channel 
flow (7.e. flow between fixed plates), Poiseuille flow, and Couette flow. It also 
includes torsional and other flows for which we do not work out the details, 
because, under reasonable body forces, such flows are not compatible with the 
dynamical equations unless inertia is neglected. 

The solutions obtained here are in terms of three unspecified real functions 
of one variable which depend on the particular material; we call them material 
functions. Our solutions are simply related. If experimental measurements are 
combined with any one of the solutions to determine the three material functions, 
then complete stress and velocity profiles can be predicted for the other flows. 
Thus, our work really establishes correspondence principles for the flows con- 
sidered. It is hoped that our results may help experimental rheologists to 
rationalize the treatment of their data by giving precise meaning to such concepts 
as ‘‘shear-dependent viscosity”’. 

If one uses a special theory, either molecular or phenomenological, to solve 
a particular flow problem of our class (say, simple shearing flow), then our three 
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material functions can be readily identified*. The solutions to the other flow 
problems can then be obtained from the formulae given here without further 
reference to the original special theory **. 

It should be pointed out that our three material functions do not determine 
the behavior of the fluid in all experimental situations, since the functional that 
predicts the response of a material to all situations may contain vastly more 
information than the three material functions introduced here. 

Some of the solutions and relations obtained here have appeared in the 
literature, but the derivations given have rested upon special assumptions about 
the nature of the fluid considered. Of these previous investigations, the most 
general is Riviin’s study of flows in fluids of the differential type [4]. Such 
fluids do not have long range memory and hence constitute only a special class 
of the general fluids considered here. It turns out, however, that RIVLIN’s 
solutions (including those for helical flow, which will be discussed in a future note) 
remain valid for general fluids. 

This paper can be read independently of reference [/] if the reader is willing 
to accept equations (1.1) —(4.3) below. 


Notation. We denote vectors and points by bold face Latin or Greek minus- 
cules: 0), 1ar, & 

Second order tensors are always denoted by light face Latin majuscules: 
A, B,I,Q. The transpose of a tensor Q is written Q*. The unit or “metric” 
tensor is denoted by J. The symbol ¢;; denotes the physical components of T 
relative to an orthogonal, but not necessarily Cartesian, coordinate system, 
except when a statement is made to the contrary. We use the summation con- 
vention throughout. The matrix of the physical components of JT is denoted 
by (TZ) or |J¢,,||, and the trace of T is written tr T. 


Partial differentiation is indicated as follows: @,=0/@y or €;=6/0x;. 


1. Basic concepts 


In a general fluid ([7], §§ 21 and 22) the stress tensor S at time ¢ depends 
on the history of the motion up to ¢ and is given by a constitutive equation 
which may be stated in the form (22.12) of [7]. In this paper we discuss only 
isochoric motions; we do not make the dependence of S on the density @ (é) 
explicit, and we write 


Silt 6 (Cp) Gis), 
Here § is a functional of the tensor-valued function C! and satisfies the relation 
Q.H(C) QF = §(2C,0") (1.2) 


for any constant orthogonal tensor Q (cf. [7], (21.7)). The function C} is the 
history of the right Cauchy-Green tensor taken with respect to the configuration 


* A summary of results pertaining to the normal stress effect for special theories 
has been presented by Markovirz [2]. The present investigation gives a proof of 
his equations (12)—(16), which he offered as conjectures. 

; ** For example, much of the work in Chapter IV of [3] could have been avoided 
if the present results had been known to the author at the time. 
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at time ¢ and is defined by 
Ci(s)=C,¢+s)=R'¢+s)R¢+s), so, (1.3) 
where /;(¢+-s) is the displacement gradient at time ¢+ s relative to the configura- 
tion at time ¢ (cf. [7], §§ 6 and 7). 
In this paper we consider homogeneous incompressible fluids. The stress is 


then determined by the history of the motion only up to a hydrostatic pressure #, 
and (1.1) must be replaced by 


T(t) = (Ci), (1.4) 
B= Spl (1.5) 


is called the extra-stress. Moreover, the tensor-valued functional is determined 
only up to an arbitrary scalar-valued functional of C}. We remove this indeter- 
minacy by the normalization 


tr tr ONC 0. (1.6) 


where 


The pressure # then reduces to the mean pressure 
=—HtrS, (1.7) 


and TJ is the deviatoric part of S. We note that the extra stress T must vanish 
if the motion reduces to a state of rest, in which case we have C;(s) =I =const. 


H 
ee. = 5iLe (1.8) 


Consider a particle which at time # finds itself at the position x, with general 
coordinates x’, and which was at time ¢+s, s<0, at the point &, with coordinates 
&'. For the dependence of & on #, t, and ¢+s, we write 


E=9,(v,¢-+s). (1.9) 


The function %,, with coordinates Oy, is the displacement function relative to 
the configuration at time ?. 

It follows from (1.3) that the covariant components of C;(s), relative to the 
natural basis at # of the coordinate system, are 


cl} «; (8) =Kyilé +3) olf +9) Ser +3), (1.10) 
hinjlt + 8) = 6; [8% (@, t+ 5)], (1.41) 


where 


and g,)(¢-+s) denotes the covariant components of the unit tensor J at §= 
%,(x, t+). If the velocity v(§, +s) is given asa function of the position § and 
the time f+ s, the displacement function 


E(s) =9,(x, ¢+ s) (1.12) 
is that solution of the differential equation 

Se ALAtt 

= v(Et+s) (1.13) 


which satisfies the initial condition 
E (0) = 9, (a, t) =a. (1.14) 


XO}: 


292 BERNARD D. COLEMAN & WALTER NOLL: 


If the motion is a steady flow, the velocity v(§) depends on only the position § 
and not on the time ¢+s. In this case 9,(a#,¢+-s), and hence also Fj(t+s) and 
C‘(s), are independent of the reference time ?. It follows from (1.4) that the extra- 
stress 7 must be independent of time. 

We assume that the body forces g (per unit mass) have a single-valued 
potential w: 


g = — grady. (1.15) 
We introduce the modified pressure: 
p=p+ oy. (1.16) 


It coincides with the ordinary pressure ~ if there are no body forces. The 
dynamical equations take the form 


div T — grady=od. (1.17) 


2. A class of steady flows 


We consider here steady flows for which the velocity field w(x) has the 
contravariant components 
fol, v2) 09} = {0, 0 (44), OF (24) 


in some appropriate orthogonal coordinate system x1, x?, x3. The matrix ||g,,|| of 
the covariant components of the unit tensor J is then diagonal. We further 
assume that the g,; are independent of x?, so that ||g;,|| has the form 


[ex(x4,23) 0 0 
le @l= 0 fexyR oY, (2.2) 
0 0 [gg (a4, x3)? 


where g,, 3, and g, are three positive functions of x1 and x3. The flows defined 
in this manner are isochoric, because the equation 


divv=0, ie, (Ve)%e,(Ver)=0, (2.3) 


where |/g = 8,283, is obviously satisfied. 


In the special case of the velocity field (2.1), integration of the differential 
equation (1.13) with the initial condition (1.14) yields the simple result 


Glee ah Cte ec u(at le 6 Fag: (2.4) 


From (2.4), (1.12) and (1.11) we obtain 


O 
IIfing + 5)|] =|] so’) 1 
O 


where v’ denotes the derivative of v. Equations (2.4) and (2.2) imply that 
fu Cras at Ce yen wl): (2.6) 


4.¢., that the covariant components of the unit tensor remain constant along the 
paths of the particles. Using equations (1.10), (2.5), (2.2), and (2.6) we obtain 
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the following expression for the covariant components of the tensor C}(s): 
[ea(2*, x*) |? + [s go(xt, x*) v'(x2)? s[go(xt, x9) Po'(xt) 0 
led I= 8 [g2 (21, 2°) ®o’(x2) [go (x4, »°)) Omen Tee) 
0 0 [ga (x*, x*) J? 
The matrix (C;(s)) of the physical components of C/(s) is 


dN (s) AE ye2ict es 
C = ef : — a v6 C 
(C:(3)) nit, x9) g(a a5 fy See UE 22) 
O O 4 
where 
xa) (2.9) 


In the remainder of this paper we shall always use the physical components of 
second order tensors. 


3. Reduction of the general constitutive equation 


The results of the previous section may be summarized in this way: If the 
velocity field has the form (2.1) in an orthogonal coordinate system obeying 
(2.2), the function C; has the form 


Ci(s) =I +sA +s?B, (3.1) 
where the matrices (A), (B) of the physical components of A and B are given by 


On AGO 4 0 0 
(A) =xl11 0 0 (B)=2||0 0 0 (3.2) 
0.0220 eh 20 


When (3.1) holds, the function C} is completely determined by A and B, and 
the functional § of (1.4) reduces to a function f) of the two tensors A and 5b, 
so that 


T=V(A, B). (3.3) * 
It follows from (4.2) that the tensor function ) must be isotropic; 1.e., 
§(Q4 Q", OBQ") = Q5(4, B) Q" (3.4) 


must hold for all orthogonal tensors Q. It is a consequence of this fact that the 
component functions §;; of are the same for all orthonormal bases. Since physical 
components are components relative to orthonormal bases, the functions },;; do 
not depend on the coordinate system or the point under consideration, if physical 
components are used. 
The normalization (1.6) of the functional implies that the function } is 
normalized by 
tothe ) 0. (3:5) 


* We note that, if (3.1) holds, A and 2B coincide with the first two Rivlin-Ericksen 
tensors A, and A, as defined in [J], (8.9). Equation (3.3) has the form of a constitutive 
equation of the differential type (cf. [7] (23.14)). 
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Furthermore, it follows from (1.8) that 


We now consider the orthogonal tensor Q whose physical components are given 
by the matrix 


1. chOey 10 
(Qo st 4 (3.7) 
0 0-1 


Then, by (3.2), we have : 
OAO =A, (ORO =.7. (3.8) 
Hence, from (3.4) we get 
QOTQ'=T (3.9) 


for T =§(A, B). Denoting the physical components of T by t;;, we find that the 
matrix of the physical components of Q TQ? is given by 


fy hie —hys 
(Q TQ") = toy too — fgg ||. (3.10) 
—ts, —tse tgs 


According to (3.9), this matrix must coincide with the matrix (7) =||¢;; , which 
is possible only if ¢,,.=t,;=hs=hs=0. Hence {/) must havesthe ton 


hi he O 
(T) =||t1 te2 0]. (3.44) 
O 0 tgs 
By (3.5), we have 
hit toe + t33 = 0. (3.12) 


From the fact that the component functions );; are independent of the 
coordinate system, and from (3.2) and (3.3), it follows that the physical components 
t;; of T are functions of x alone. It is clear from (3.11) and (3.12) that these 
components ¢;; may be expressed in terms of three independent functions of z. 
We choose 


ho=t(x), by—tg=0,(%), bey —ty3 = 02(x). (3.13) 

For the physical components s,;;=1;,—0,; of the stress tensor we also have 
Sig = Soi = T(x), (3.14) 
S11 — $33 =01(%), — Sg3 — S33 = 0g (x). (3.15) 


The three functions t, o,, and o, depend only on the material, not on the particular 
flow; they are material functions. 


The functions 7, o,, and o, are not completely arbitrary, as may be seen from 
the following argument: For a second choice of Q in (3.4) we take 


10 0) 0 


ee es NG) (3.16) 
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Equation (3.2) then yields 


OAO = Ae OBO = 8, (3.17) 
and, from (3.11), we get 
Z ty —he 0 
(were = —ty, tye 0 (3.18) 
oat 6 


It is clear from (3.2) that the conjugation operation (3.17) is equivalent to 
replacing x by —x. Thus, by (3.13), equation (3.18) states that 


t(—m%) =—t(x),  9;(x) =9;(—%); (3.19) 
1.e., that t must be an odd function, while o, and o, must be even functions. No 
further restrictions are imposed on 1, o,, and a, by the isotropy condition (3.4). 


The condition (3.6) implies that the material functions 1, o,, and o, must 


vanish for x =0: 
T{0)) =, (0) = cH(O\i= OF (3.20) 


A further restrictive condition on t follows from the principle that the 
dissipation of energy must be positive. This condition is 


eT ON uae = (> (3.21) 


1.e. T(x) must have the same sign as x. If it is assumed that t is twice con- 
tinuously differentiable, then (3.20) and (3.21) imply that t(x) must be a strictly 
increasing function of x in some interval —x,S*xS+%, around x=0. In this 
interval t will have a strictly increasing and odd inverse, tT. In the remainder 
of this paper, whenever T occurs, we assume that we are in the range in which 
T is defined. If 7 is not defined over the entire real axis then the flows which 
we discuss below will not occur at high speeds. 

We can give a precise meaning to the “‘shear-dependent viscosity’ used in 
the rheological literature by identifying it with the material function 7: 


Fie) (3.22) 


Clearly 7 must be an even function. Since t(0) =O and 7 is twice differentiable, 
y is differentiable, even at x =O, and 


n'(0) = 0. (3.23) 
If we assume that o, and o, are differentiable, we also have 
0; (0) =03(0) =0. (3.24) 


If a polynomial approximation for y, o,, and o,, in a neighborhood of x =0, is 
used to fit experimental data, only even powers of % can occur. 


4. Rectilinear flow 


If the velocity field has the form (2.1) in a Cartesian coordinate system 


(i= Cae eS Wnte hoes (4.1) 
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then we call the motion rectilinear flow. In this case 


Si Ga ese (4.2) 
and hence, by (2.9), 
HU (4.3) 
By (3.13), the extra stresses ¢,, are then functions of x only, and the dynamical 
equations (1.17) become 


Orter —49=0, O%t,—-A97=0, O,9=0. (4.4) 

A simple analysis shows that these equations are satisfied only if 
iO) = — GeO, (4.5) 
P=9(%,9) =—ay +4, ,(%) —¢, (4.6) 


where a, b, and ¢ are constants. 


There are two important special cases of rectilinear flow which we shall 
discuss in the next two sections. 


5. Simple shearing flow 


This is a rectilinear flow between an infinite plate I at rest and an infinite 
plate II moving with a constant speed V parallel to plate I. We denote the 
distance between the plates by d. We choose the y and z axes to be in the plane 
of plate I and such that plate II moves in the y direction. The % axis is then 
perpendicular to the plates with x =0 at plate I. Assuming that the fluid adheres 
to the walls, we have the boundary conditions 


DO 0; bed ae (5.4) 


We assume now that the modified pressure g, defined by (1.16), has no 
gradient in direction of the flow; 7.e., that 0,q@—=0. We then have a =0 in (4.6). 
Hence, by (3.13), (4.3) and (4.5), the shearing stress 


by y (X) =e (VX) =O (5.2) 


is constant. It follows that v’(x) = (6) must be constant and hence that v(x) 
is linear. The conditions (5.1) then imply that 


V 
v(t) = 7G (5.3) 
and (4.3) becomes 
V 
“=~, = const. (5.4) 


According to (3.14), the shearing stress is 


veeilna ee 


Since Spy is equal to the tangential force per unit area that must be applied to 
the moving plate IT in order to produce the flow, equation (De 5) can be used to 
determine the material function 7 experimentally. 
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It follows from (5.4) and (3.13) that the ¢,, must actually be constant. Hence, 
by (4.6), since a =0, 
Grate 6 == CONS: (5.6) 


Using (1.16) and s,;=t,,—, we find 


and (3.15) becomes 7 
Sipess, hay (=): Syy — S15 92( 3). (5.8) 


Hence, the matrix of the stress tensor is 


c— oy (3) 0 | 


6) O c—op—o,(~ 


Since the normal stresses can be measured, at least in principle, equations (5.8) 
can be used to determine the material functions o, and og. 


6. Flow through a channel 


This is a rectilinear flow between two parallel infinite plates which are both 
at rest. We now place the y and z axes in a plane half-way between the plates; 
the x axis is then perpendicular to the plates. Let the distance between the 
plates be d; it follows that the plates are in the planes x = +$ d. 


The fluid adheres to the plates. For rectilinear flow, to get any motion at 
all we must have a =+ 0 in (4.6). When a is positive, there is a flow in the positive 


y-direction. From (3.13), (4.3) and (4.5) we get v'(x) =t (—ax+0). The 
boundary conditions are v(+3d)=0. The fact that t is an odd function 


implies, 0 =O; 2.¢:, a 
CAN TeX). (6.1) 
After integration we obtain the velocity profile 
v(x) = — ft (aédé. (6.2) 
—hd 


The total applied force f in the direction of the flow, exerted on a column 
of fluid with cross-section .7 and lying between two planes y=y, and yv=y,,, 
is given by 


yn 
be i SyyG@ A), oe (J 2) yd A), a P (2,y) odA dy. (6.3) 
This reduces to 
£= J {(Syy — OV) y= — (Syy — OP)y=u} G4 - (6.4) 
A 
By (1.16) and (4.6) we have 


Sy) — CY =byy — P — OY = by — 9 = 2 + ty (4) —4,, (4) Pe, 
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hence Syy — Op =ay+h(x). (6.5) 


Substitution into (6.4) yields 
f=a(vy—M)A, 
where A is the area of the cross-section of the column. We thus have 
if 
i 6.6 
A(Yq1— 1) ’ e0) 
which gives a physical meaning to the parameter a. It is the applied force per 
unit volume in the direction of the flow. It is the driving force that produces the 
flow. If there are no body forces, a reduces to what is usually called the pressure 
head per unit length. It must be produced by a pump or similar device. If the 
only body force is gravity, acting in the direction of flow, and if there is no 
pressure head, then a is simply the specific weight of the fluid. In any case, 
a is a measurable quantity. We call it the driving force. 
The volume discharge per unit time through a cross-section of the channel 
one unit deep is 
4d 
Q= Jf v(x)dx. (6.7) 
—id 
Integrating by parts, remembering the boundary conditions v(+ $d) =0, and 
using (6.1), we find 


Q=— f xv'(adn= f xtlaxdz. (6.8) 
ia —thd 
We finally obtain 
sad 
Q=— | Ez) aé. (6.9) 


Thus, the experimentally measurable rate of discharge Q is related to the 
measurable driving force a through the material function t which can be deter- 
mined by independent measurements. 


Equation (6.9) is easily solved for the function 7 in terms of the function 


Q=Q(a): 


-l/ad Dy Gl 
7 S\= oa <—[a? Q(a)]. (6.10) 
This formula may be used to calculate the material function t when Q has been 
measured as a function of a. 
The normal stresses are obtained from (4.6), (3.15), (4.3) and (6.1): 


Sex =4Y + QY—E, 
Sye— Sz, =0,(t (a2), (6.11) 
Syy> Sey == Oe (z (a x). 
If the material functions 1, o,, and o, are known from independent measurements, 
these formulae may be used to predict the measurable normal stresses. It will 


be noticed that the normal stresses are determined only up to a constant hydro- 
static pressure of magnitude c. 


a 
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If only t is known, then (6.11) may be used to determine o, and o, by meas- 
uring s,,—s,, and s,,—s,,. If the potential y of the body forces is known, then 
a measurement of s,, just determines the value of the constant c. 


7. Poiseuille flow 


We consider a steady flow through an infinite circular pipe of radius R and 
use cylindrical coordinates with the z-axis coincident with the axis of the pipe. 
We label the coordinates as follows: 


PEE BES SSE: (7.1) 


We say that we have Poiseuille flow if in this coordinate system the velocity 
field has the form (2.1) and if the fluid adheres to the walls of the pipe; 7.e. 


v(R) =0. (7.2) 


The coordinate system 7, z, # is orthogonal and obeys (2.2) with 


Si ba es = (7.3) 
Therefore, we find that (2.9) becomes 
fata (7.4) 


By (3.13), the extra stresses ¢,; are then functions of 7 only, and the dynamical 
equations (1.17) become 
Or a tas (t,, me ty) =< a, =0, 
0b ert =e p =); (7.5) 


— Op =0. 


An elementary analysis shows that these equations are satisfied only if 


t,,= —gar+obri, (7.6) 
P= Plr,2)=—az+h(7), (7.7) 
h'(r) = Ont, i us (t,, 3 tye) , (7.8) 


where a and 0 are constants and /1 is an arbitrary function of 7. Since the shearing 
stress s,,—#,, must be continuous at the center of the pipe (7.e. for y=0), we 
must put b=0 in (7.6). According to (3.14) and (7.4), we thus have 


hence, since t is an odd function, 
=] - 
V(t (5 a7) (7.9) 


This equation, together with the boundary condition (7.2), determines the 
velocity profile: 


v(r) =f T(bad) dé. (7.10) 
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An investigation analogous to that given in § 6 for the case of a flow through 
a channel shows that the parameter a is again the applied force per unit volume 
in the direction of the flow. We again call a the driving force. It reduces to the 
pressure head per unit length if there are no body forces. If the flow takes place 
in a vertical pipe under the action of gravity alone, a is just the specific weight 
of the fluid. 

The volume discharge per unit time through a cross-section of the pipe is 


R 
Q=2x f v(r)rdr. (rs) 
0 
After integration by parts and use of (7.2) and (7.9), we find* 
4aR 
_ 8n f Ps 42 
= as plc AC (7A2) 
0 


Thus, if the material function t is known, one can predict, for pipe flow, the 
dependence of the rate of discharge Q on the driving force a. 


As in the case of channel flow, (7.12) is easily inverted: 


@ (aR) = aps Gg (O(a). (7.13) 


If Q=Q(a) has been measured, this formula may be used to calculate the 
material function T. ; 
To calculate the normal stresses, we observe that (3.13), (7.4), (7.9) and the 
fact that o, is an even function imply 
~~ tea = (7 Bar), (7.14) 
Integration of (7.8) yields 


R 
h(r) =1,,— f €40,(7 $aé)) dé —c, (7.15) 


where c is a constant. Substituting (7.15) into (7.7) and using p=p+oy and 
S,,=1t,,—P, we obtain 


R 
Sy = 42+ oy +f &40,(z aé)) dé +c. (7.16) 
By (3.15), we have 


$7 — 899 =0,(t (Aan), Sez — S99 = 0n(t (a7). (7-47) 


The remarks made in §6 about the prediction of measurable normal stresses 
in channel flow also apply here. 


8. Couette flow 
Here we consider a steady flow between two infinite coaxial cylinders; the 
inner cylinder has radius R,, and the outer cylinder has radius R,. Using cylin- 
drical coordinates with the z axis coincident with the axis of the cylinders, we 


* Equivalent relations have appeared previously in the literature (cf. [5], pp. 249 
and 250). The derivation given here, we believe, is the first that is not based on special 
assumptions. 
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now label the coordinates 
Pie Pa fee 4, (S24) 


If, in these coordinates, the velocity field has the form (2.1), we speak about 
Couette flow. We than have 


so that (2.9) becomes 


Note that v(7) has now the physical dimension of an angular velocity. By (3.13) 
and (8.3), the extra stresses ¢;, depend only on 7. The dynamical equations (1.17) 
become 


One zie ip (t,, a ty») a. oC, ar C12 
0,t,9 + 2741, — 1 Bg =0, (8.4) 
Oi Oe 


An analysis of these equations, using the fact that m must be single valued and 
hence periodic in #, shows that m can depend only on 7 and that t,,=M/(227°), 
where M is a constant. From (3.13) and (8.3), it follows that 


La= tev") =—— (8.5) 


The moment per unit height exerted on the fluid inside the cylindrical surface 
vy =const. is, by (8.5), equal to 
r (2207) t,y=M (8.6) 


and is independent of 7. From (8.5) we get 


vin=—t| 


pacslen ; (8.7) 


227 
We consider further the special case in which the inner cylinder remains at 


rest while the outer cylinder rotates with constant angular velocity 2. Assuming 
that the fluid adheres to the walls, we then have 
OG) — Om OUR) = 2. (8.8) 
M is the moment per unit height that must be applied to the outer cylinder 
in order to produce the flow. 
Integration of (8.7) with the boundary condition (8.8) yields 
M/(2 x R?) 
— 1 we 
eS | 
M/(2 x R2) 


<1 (6) dé. (8.9) ** 


Thus, a knowledge of the material function t makes it possible to predict the 
angular velocity 2 of the outer cylinder as a function of the moment per unit 
height M. If the gap between the two cylinders is very small [(R,—R,)/Ri< 1], 


* Note that this differs from the labeling (7.1) used for Poiseuille flow. 
x* A derivation, based on special assumptions, of a formula equivalent to (8.9) is 
given in [6], § 3. 
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the formula (8.9) may be approximated by 


Re— he et M ' 
Lad — a 
Dw Lt (ee) (8.10) 


We wish to invert (8.9) so as to be able to calculate the material function t 
from an experimental determination of 2 =Q(M). Differentiation of (8.9) with 


respect to M gives 
dO), ete a M ae 
esa ee (sani) Desa tee) 


With the abbreviations 
QM) _ Ri 


F(M) =2M 2 ar RE (8.12) 
equation (8.11) becomes 
= Me ol eee 
SACU Aes a (8.13) 
If we replace here M by x” M and sum from ~=0 to n=N, we get 
af o” es M Sif estes | TULL 
2th M7 sear Tt (> ae at ‘ (8.14) 


Since R,<R,, we have «<1 and hence lim «&”*1—0. The continuity of t at 
n—> co 
zero and the fact that 7 (0) =0 (cf. (3.20)) insure that 


Sah fa M 
ingle Se 
Hence, (8.14) yields* 
-1/ M Shae 
T faa = > Fe M). (8.15) 


Equations (8.12) and (8.15) make it possible to calculate t if 2 =(2(M) is known. 
When the gap between the cylinders is small, t can be calculated directly from 
(8.10). 


We now determine the normal stresses. From (3.13), (8.3), and (8.7), we get 


br toamen(? (ga) ea (aa) 19 


Substitution of (8.16) into the first equation of (8.4) and integration yields, after 
using SP and Y =p oy, 


sen fete ba Ma) —Le aM ae res 7 


1 


where v(&), by (8.7) and (8.8), is 


v(é) = ue T ( sam aC. (8.18) 


* Note added in proof: It has been brought to our attention that an equivalent 
inversion of (8.9) was published by Krimecrer & ELrop [7]. 
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By (3.15), we have 


-f lp IGE =y 
Sy1— Sy, = 04 (7 (=) , S00 8.2 =02(7 | 


saat))> (8.19) 


The remarks made at the end of § 6 apply here also. 


When the potential y of the body forces is independent of 7, (8.17) yields the 
following expression for the difference of the normal stresses at the outer and 
inner cylinders: 


Ry 
te sexbht Aske ValieNeN Went dy o.(-lt aie ™ 
S; (Re) —s, (Ry) = ( {o r{v (7) |? a yA a,(z ian ery a(t (saa) har (8.20) 
Rk, 
When (R,— R,)/R\<1, this formula may be approximated by 
ee Py cated Eye 
S, (Ra) — 8, (Ra) & HAE lo. (z eae ) o,(z Cea ))}- (8.21) 


The error in both the approximate formulae (8.10) and (8.21) is of the order 
of [(R,—R,)/R,]?.. The inertial term o7[v(r) |? in the integral on the right side 
of (8.20) makes a contribution to s,,(R,) —s,,(R,) of the order of [(R,— R,)/R,]* 
and hence does not appear in (8.21). 
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The Constitutive Equations 
for Certain Classes of Deformations 


RSS RVLIs 


1. Introduction 


In previous papers [/, 2] constitutive equations have been developed for a 
certain class of deformations in materials for which the stress at any instant 
of time depends, in general, not only on the displacement gradients at that 
instant of time, but also on the previous deformation history of the material. 

We considered [/] that a material is undeformed until time zero and that, 
in the short time interval 0 to ¢,, it is subjected to a deformation which is there- 
after maintained constant. It was assumed that the material is such that for 
deformations carried out in a short time ¢,, the material behaves substantially 
as a perfectly-elastic material, but for deformations carried out during much 
longer times the stress depends on the deformation history up to the instant 
considered. For the limited class of deformations which are maintained constant 
after time ¢,, it is apparent that the stress at an instant of time ¢(>>4,) will depend 
only on the (constant) displacement gradients and on the time ¢. With this 
as starting point, it was shown that constitutive equations closely related to 
those which are valid in finite elasticity theory can be obtained. 

For example, in the case when the material is incompressible and isotropic 
in its undeformed state, the constitutive equation takes the form 


i p 0;; Oya Ie Coe eee ea (1.1) 
Here g;; is defined by 
_ Ox; Ox; 
bi GES OR em 


x, and X; are the coordinates in a rectangular Cartesian coordinate system x 
of a generic particle of the body in the (constant) deformed state and in the 
undeformed state respectively and ¢;; are the stress components in the coordinate 


system x. In equation (1.1), £ is an arbitrary hydrostatic pressure and a, and a, 
are functions of ¢ and of J, and J, defined by 


I,=g;; and To = (8:18); — 815 8a) - (1.3) 


It was then shown that formulae can be derived for the forces necessary to 
maintain certain simple physically important types of deformation of the class 


considered, without making further assumptions regarding the manner in which 
a, and «, depend on J,, 7, and #. 
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In a subsequent paper [2] it was pointed out that the assumption that the 
material behaves as a perfectly-elastic material for deformations carried out 
during short times ¢, is an unnecessarily restrictive one for the validity of equation 
(1.1) as the constitutive equation for deformations of the class under consider- 
ation in an incompressible, isotropic material. We may instead merely assume 
that the stress at time ¢ has a ‘“‘memory”’ of the deformation to which the material 
was subjected at earlier times which decays sufficiently rapidly. 

It was also shown [2] that if the gradients of the (constant) displacement wu; 
are small compared with unity, (1.1) may be re-written as 


tj = — 2055+ Bi es; + Be Cin en; (1.4) 
where e;; is defined by 
welahp fits Ou; (1.5) 
N= 2 lee r we : 
and f, and f, are functions of ¢ and of J, and J, defined by 
Fi = 3 (655% 1 — %158%;) and J,=|e;,]. (1.6) 


Furthermore, the validity of such equations for a certain class of materials 
was demonstrated. It was found that for these materials 6,=0 and f, took 


eke EE UAR (1.7) 


In the present paper, we shall show that constitutive equations which are 
formally similar to (1.4) may be derived on the basis of the assumption that 
the displacement components in the coordinate system »% all vary with time 
according to a specified law which is the same for each of the displacement 
components and independent of the particle considered (see equation (2.3) 
below). The forces which must be applied to the ends of a rod of circular cross- 
section in order to produce in it simultaneously a torsion and simple extension 
for which the amount of torsion and fractional extension vary with time in 
the same manner are then calculated. Finally, we discuss the case when the time 
dependence of the constitutive equation occurs as a common multiplier of all 
the terms in it. 


2. Basic Assumptions 


We consider a body to undergo a deformation such that a particle of the 
body which is initially (at time t=0) at X; in the rectangular Cartesian coordinate 
system x moves at time t to x, in the same coordinate system. The displacement 


components w; in the coordinate system x are defined by 
U;,=%,—X; (2-4) 


and the deformation is completely specified if the functional dependence of wu; 
on X; and 7 is given by a relation of the type 


u;—=u;(X;,T). (2.2) 
Now, let us consider the class of deformations for which the relation (2.2) has 
the form Finan 233) 


where /(t) is some specified function of t and v; is a function of X; only. 
Arch, Rational Mech. Anal., Vol. 3 W 
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We shall now assume that once we have specified the function f(t), the 
components of stress ¢;; at time ¢ in the system x at a point of the body, which 
was at X, at time t=0, depend only on the gradients 0v,/0X; at the point; 7.e. 


t;; = Pi; (Ov,/OX,, t), (2.4) 


where ;, is a single-valued function of the arguments 0v,/0X, and ¢, which 
is completely determined by /(r). 

We now assume that @u,/0X,<1. It is then possible to impose certain 
restrictions on the manner in which g;; can depend on dv,/0X,. These are 
discussed in Section 3. 


3. Derivation of the Constitutive Equation 
Let us now consider that the body undergoes, in addition to the displacements 
defined by (2.3), small rigid rotations of amounts /(t){,; about the axes ~%,. 
Let u; be the displacement components of a generic particle of the body, initially 
at x;, in these rotations. Then, provided the rotations are sufficiently small 
so that 0u;/0x;<1, we have 
U,= f(t) &5 542; Xp- (3.1) 


Since we have already assumed 0u;/0X ;<1, it follows, with (2.1) ,that 0u,;/0X;<1. 
Neglecting terms of the second degree in du,/0X; and 0u,/0X; in comparison 
with those of the first degree, and employing the notation 


0; = €;;,.2; X,, (3.2) 
we have 
(2) = i Ds 
aoe HN eet tage (3.3) 


With the same approximation, we may take the stress components at time f 
in the coordinate system x associated with the displacement field u,;+7#; to be 
the same those as associated with the displacement field u;. We thus have, 
with (2.4), 


OUp ise OUp OV, 
vil aX,’ j vil ax, v Ox, i) 6-4) 
for alld, of the form (3.2), 7.e. for all 0v,/0X, of the form 
Dp 
gare oeagedee (3:59 
Taking 
1 C0, 
Q,= 2 Eymn ry or, (3 6) 
we obtain, from (3.5), 
OUp 1 OUm 
ax iy oF ®rmn rpg ox Od) 
Since 
Eymn Evpg = Om p Ong One Ons (3.8) 


equation (3.7) yields 


Certain Constitutive Equations 307 


From (3.9) and (3.4), we obtain 


OUp oi 1 / OUp Og 
vil ax. ) vil (ax, cc aX)’ i). oe 
Since 
aU, Ob, 
| Gs 
Dy er eee (3.11) 


for all dv,/0X, of the form (3.5), the validity of the relation (3.10) implies the 
validity of the relation (3.4) for all 0v,/0X, of the form (3.5). We thus have 
the result, from (2.4), that 

ti7= Vij lpg. 4), (3.12) 
where ¢,, is defined by 


ii OUp Og 
Gar lan: a Bree 3.13) 


4. The Isotropic Case 


We now consider that the material of the body is isotropic at time t=O. 
Let x’ be a generic rectangular Cartesian coordinate system related to the system x 


b ’ 
where a;,; are constants satisfying the orthogonality conditions 
Ajj; Az, = 4; 5 Ay, = 0,5. (4.2) 


Let X; be the coordinates in the system x’ of the point which is at X; in the 
system x. Let ¢;; and u; be the stress components and displacement components 
respectively in the system x’. Then, we have 


, / / 5 
Xi =AjypXp, Uj =AiyUy and tb; = ajy aja ty. (4.3) 


Employing the notation ui, = f(t) v;, 


we obtain from (2.3) and the second of equations (4.3) 


v;= Bip Yp- (4.4) 
Defining ¢,, by os et 
r __ 4 (OU, “1 2 
OAS ae oh ax]? 2) 
we see that 
Ci; = Vip iq pq: (4.6) 


It follows from (3.12) and the fact that the material is isotropic at time 
= th t , , 
ares ah tis = Pij (Cpa 4): (4.7) 
From (3.12), (4.3) and (4.7), we obtain 
Aim in Pmn (pqs t) = Piz (Cbg 4) (4.8) 
where ¢,,, is given by (4.6), for all a,;; satisfying (4.2). This equation represents 
the limitation imposed on the form of ,; by isotropy of the material when 
T=0. 
21* 
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Let y,; and y;; be the components in the coordinate systems x and % re- 
spectively of an arbitrary symmetric Cartesian tensor of the second order. Then, 


Vann — Aim Fin Piz: (4.9) 


Multiplying the relation (4.8) throughout by y;,;, and employing the relation 
(4.9), we obtain 


Pii Pij (pq 4) = Wij Pij(Cpq ¢)- (4.10) 
With the notation F=y,; 9; ;(@pq, ¢), we see that 
Nf Ole oF 

ny : 4.11 

YOO p oe a) ( ) 


The relation (4.10) expresses the scalar invariance of F under orthogonal 
transformations. 

If we assume initially in (2.4) that @,;,; is a polynomial in 6v,/0X,, then in 
(4.10) p;,; is a polynomial in e,,. F is then a polynomial in e,, and y,,, linear 
in the latter. / must therefore be expressible as a polynomial in the elements 
of an integrity basis under the orthogonal group for the two tensors w,;; and e;;. 
With the notation 

b=|lyii|| and e=|les,ll, (4.12) 
we thus obtain 


F=a,trp+a,trpe+a, trope’, (4.13) 


where a, «, and a, are polynomials in tre, tre? and tre® and are single-valued 
functions of ¢. From (4.13), (4.11) and (4.12), we readily obtain 


bij = Aq Of; + Oy C5; + Oy Cry Cp; (4.14) 
If p;; is a single-valued function of e,,, then we can employ a procedure 
used, in a somewhat different context, by EricksEN & RIVLIN [3] to show 


that an equation of the form (4.14) is still valid with a», «, and a, single-valued 
functions of tre, tre”, tre’ and ¢. 


5. The Incompressible Case 
If the material considered is incompressible, we have, still assuming @,/ OX ;<1, 
Ou; 


ae Os (5.4) 


With (2.3) and (3.13), this equation may be rewritten as 
tre=0. (5.2) 


Also, noting that if the deformation is specified the stress is undetermined to 
the extent of an arbitrary hydrostatic pressure — p (say), we replace our initial 
assumption (2.4) by the relation 


Ox, 
tis Pril axe i P 95;- 63) 


We then obtain in a manner similar to that adopted in Sections 3 and 4 that, 


in general, 
bij = Piz pq, t) — 9 0,; (5.4) 
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and, if the material is isotropic at time t=O, 
bj = — £05; +01 C55 + Oy C jp ep ;, (5.5) 


where « and a, are single-valued functions of tre?, tre® and ¢, and are poly- 
nomials in tr e? and tr e? if y;, in (5.3) and (5.4) are polynomials in the arguments 
0xp/0X,. 


6. Derivation from the Constitutive Equations 
for Materials with Memory 


We consider a body to undergo the deformation described by (2.2). We 
assume that the material of which the body consists is such that the stress 
components ¢;; at time ¢ are single-valued functions of the displacement gradients 
at time ¢ and continuous functionals of the displacement gradients in the time 
interval t=O to t=#. Then, the stress components can be expressed with any 
desired accuracy by an expression of the form 


fh? + fa (t, T) ) ate dt, + 


Ou 0 
~ / ip Pasir (ts Th, Te) “URI Shel dey day + (6.1) 
oo ODD. 
where i ie jfl,»--+) 4n addition to being continuous single-valued functions of 


the indicated variables are also single-valued functions of du, (t)/dX,. 
Now, if we consider only the class of deformations described by (2.3), with 
specified /(t), we see that equation (6.1) becomes 


a) 
= M+ oe “Ah (Pan (tT) fe) dea + 


one f / Wasnt (ts Tas Te) (ta) (ta) dtd t% + (6.2) 


Stats 


Carrying out the integrations, we see that ¢;; may be expressed as a single- 
valued function of dv,/0X, and ¢, as indicated in (2.4). 

From similar considerations it is evident that if we take as our starting 
point that ¢;; are single-valued functions of the displacement gradients and 
their time derivatives of the first R orders (say) at time ¢ and continuous function- 
als of these in the time interval t=O to ¢, we arrive at a similar conclusion. 


7. The Solution of a Simple Problem 


If we have a body of the material considered and we know the constitutive 
equation (4.14) or (5.5), for a given function /(¢) in equation (2.3), then we can 
always calculate the forces necessary to produce in the body a specified defor- 
mation of the form (2.3), without further assumption regarding the precise 
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manner in which the «’s depend on the J’s. In general, however, the calculated 
force system will involve body forces which could not, in practice, be readily 
applied, or will lack physical interest in some other respect. 

However, in certain simple cases it can be shown that the deformation may 
be produced by surface tractions only and these may be calculated. As an 
example, we may consider the forces necessary to produce, in a rod of circular 
cross-section with radius a of incompressible material for which the constitutive 
equation for the class of deformations (2.3) has the form (5.5), a simultaneous 
small extension and torsion. 

We take as our reference system a rectangular Cartesian coordinate system x, 
the x axis of which coincides with the axis of the rod. If the fractional extension 
at time 7 is ef (t) and the amount of torsion is p/(t) and ef(t)<1 and y%3f(t)<1, 
we may write the displacement components w; in the rod at time ¢ in the form 
(2.3) where v, is given by 

Vv, = —pX,X,—FeEX,, vo =ypX,Xz—FEeXs, Vg=—EXsz. (7.1) 


We note that the relation (5.1), which expresses the fact that the volume of 
each element of the body is unchanged in the deformation, is satisfied. Sub- 
stituting (7.1) in (3.13) and the equations so obtained in the constitutive equation 
(5.5), we see that the stress components at time ¢ are given by 


bi=—p-— FH e4 ZX (e2-+ y? X92), tog = 2 YA1(% + 3 EM), 
tog = —p—FHe4 ty (€? + yp? X4), tay = —ZpXeq(% + 3 EM), (7.2) 
tyg = — pt oy e+ a (6? + 4 R?) he= — Zp XX. %, 


where 


RSX -P x 
In these expressions, «, and a are, of course, functions of ¢ and of tre? and 
tr e® given, from (7.1) and (3.13), by 
tre=se+5y?R? and tre’=2e3+ sey? R2. (7.3) 


If the function /(t) varies sufficiently slowly with time, the deformation 
may be treated as a quasi-static deformation and inertial forces may be neglected 
in the equations of motion. With the approximation indicated by 0u,/0X,<1, 
we may then write these as 

Otis 

OX; = O. (7.4) 
With (7.2), and the assumption that the surface tractions are zero on the cylind- 
rical surface of the rod, equation (7.4) yields 


R 
p= — (oe Loge) Ly? fRAR. (7.5) 


Introducing this expression for p into (7.2) we find that the surface tractions 
per unit area which must be applied to the ends of the rod in order to maintain 
the deformation have azimuthal components 


2 VR (oy + 3 € ae) (7.0) 
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and longitudinal components 


R 
Soe + faze + hoppy Rt iy? f aRAR. (7.7) 


The distribution of surface tractions (7.6) and (7.7) over an end of the rod 
is statically equivalent to a couple of magnitude 


my { R8(a + $99) dR (7.8) 


and a tensile force of magnitude 


mf (30 E+ 3 a +40 y?R2) RAR. (7.9) 
0 


8. A Particular Form of the Constitutive Equation 


Let us suppose that the form of the function f(t) in (2.3) is specified, and 
that the constitutive equation (4.14) [or, in the particular case of an incom- 
pressible material, (5.5)| is valid for this form of f(t). Now, let us specify a 
deformation of the type (2.3). In general, if a certain force system is appro- 
priate for the maintenance of the deformation at one instant of time, a force 
system of quite a different type must be applied at other times. For example, 
if at one instant of time the deformation may be supported without application 
of body forces, then this will not generally be the case at other times. 


However, if %), «, and « have the forms 


=P) Yo, UM=Plt)y, and «= P(t) ye, (8.1) 


where yy), y, and y, are independent of ¢, then, if the deformation can be main- 
tained at one instant of time without application of body forces and with boundary 
conditions which take the form of specified displacement components or zero 
components of surface traction on portions of the surface which are specified, it 
can easily be seen that this will also be the case at all other times, provided tha 
that the deformation may be treated as quasi-static. For, the equations of motion 
and the boundary conditions may then be written in a form which is independent 
of ¢. 
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On Circulation in Relativistic Hydrodynamics 


A. H. TAuB 


Communicated by G. C. MCVITTIE 


Abstract. In this paper it is shown that the relativistic analogue of the circulation 
in classical hydrodynamics is the integral 


c= (1 +4) U,d x%, 


where i is the specific enthalpy of the fluid, U, are the covariant components of the 
four velocity of the fluid and the integral is taken around a closed curve traveling 
with the fluid. This integral is similar to an integral given by LICHNEROWICzZ, but differs 
in the factor multiplying U, in the integrand. 


The necessary and sufficient conditions for C to vanish are shown to be that 
the flow be irrotational and isentropic. The necessary and sufficient condition for 
C to be a constant of the motion is that the pressure be a function of the density 
alone at every point in space-time occupied by the fluid. This condition is shown 
to be violated when a shock is present. Results for similar integrals taken around 
vortex lines are also obtained. The relations between these results and Bernoulli’s 
theorems are discussed. 


1. Introduction 

It is the purpose of this paper to discuss an integral formed from quantities 
describing the relativistic behavior of a perfect fluid. This integral is analogous 
to the integral defining circulation in classical hydrodynamics. The results 
obtained below are, in part, similar to results given by SYNGE [J] and reported 
by LicHNEROWwIcz [2] but are obtained by using a different formulation for 
the stress-energy tensor of a perfect fluid. Our results hold for both special 
and general relativity and depend only on the assumed nature of the stress- 
energy tensor and the four conservation laws it obeys and the equation of con- 
servation of matter. 


A perfect fluid is described by the stress energy tensor 


Ter 6 OF Pine aee (1.1) 
where 
UU, =1; (4.2) 


U* represents the four dimensional velocity vector of the fluid, 8,», the metric 
tensor of space-time, #, the pressure, and 


ctmelttet Fr 
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In this equation @ is the proper density, and 


e= (Pp, 0) (1.4) 


is the proper, specific internal energy of the fluid. Equation (1.4) is sometimes 
referred to as the caloric equation of state of the fluid. 

The equations governing the motion of the fluid are the formulation of five 
conservation laws, namely, the conservation of mass, energy and momentum. 
Where there are no singularities, these equations take the form of partial differ- 
ential equations which are 

i os =, 
(@U"),,=0, 
where the semicolon denotes the covariant derivative with respect to the metric 


tensor g,,. This paper is not concerned with the determination of this tensor. 
It is assumed to be known. 


(1.5) 


Across a singular hyper-surface [3], the conservation laws are 
[7,4] =0, 


6 
[o U¥ 1,) =0, ee 


where 4, are the components of the unit normal vector to the singular hyper- 
surface, and we have used the notation 


[/] = lim [f(x" — 06") — f(x" + e&)J=h—h, 
where the x“ are the coordinates of a point on the singular surface and the &“ 
are arbitrary. 

If equations (1.1) are substituted into the first four of equations (1.5), they 
may be written as (cf. [3]) 

7 epee Pav 1 a v 
Ct P= a (g GEO") (1.7) 

and 


TS.,,U" =0, (1.8) 


where T is the proper temperature, and S is the specific proper entropy. The 
quantities T and S are defined as in classical hydrodynamics by the requirement 
that 


TdS =de+pa(_). (1.9) 
We define the specific enthalpy as 
(ee = (1.10) 
Th 
ey die pee eee (1.11) 
ag c= 9(c2 +7) (1.12) 
We further define b 
dp 
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where the integrand is considered a function of p and S, and the integration is 
carried out with S=constant. Then 


P 
Katt elt Senet Uda 
= o= | a OM aS GI, Sie ee 


We shall treat p, as independent of the coordinates of points of space-time in 
our subsequent discussion. Then 
Og ep on OS 
- Oye gt Ba ae ews 


(1.14) 


where p is kept constant in the partial differentiation of g with respect to S. 
For some functions ¢(p, 9) there exists a fy such that 


Pes See Onn 
With this value of py, equation (1.13) becomes 


vs 1 
é ie es (1.15) 
ge LOY Sete (ROOD ead, 
as OSes ame2 


2. The Tensor Q,, 
We define the antisymmetric tensor 
QV, c= (2.4) 
ts Vier O (2.2) 


and @ is given by equation (1.13). This tensor will play an important role in 
our subsequent discussion. We, therefore, discuss some of its properties. It 
follows from equation (2.1) and (2.2) that 


Loe =e? ha Aes cy (P., U, SP U,) | , 


where 

Wee Urea Ug (2.3) 
and hence 

Wy U” a ae Cie (2.4) 


in view of equation (1.2). 


We further define the vorticity vector v“ by the equation 


Uf eae ee U,, (2.5) 
where 
LvOT 1 vOT 
| aie a s# (2.6) 
and 
EL vce= Vee 2 Sxsars (2.7) 


The relative tensors ¢,,,,=e“"’* are zero unless all indices are different. In 
that case, they are plus or minus one, depending on whether the values taken 
by the indices are an even or an odd permutation of 1 2 3 4. 
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It is a consequence of the definition of these quantities that 


Weer LER Ag ser Apiee Ay...4 
E Ae 2; re NEAL 2 =i Al ae (2.8) 


ine DIV ras 


where the tensor One is zero unless the indices /,... A, are an even or an odd 
permutation of the indices »,...»,. In the former case it has the value +41, 
and in the latter case, the value —1. 


Multiplying equation (2.5) by E and summing, we obtain 


uaBy 
WE apy = 2 Onpy Wor U, = U, Wey + U, Wa + U, Wy 9: (2.9) 
It then follows that 
Wypm= 08 GW See on OCR UU Use, Oss (2.10) 


Substituting equation (1.7) into equation (2.10) and the resulting equation 


into the expression for 2, in terms of W,,, we obtain 


Ore #0 OU Faye (5, Ue 5,5 U,) (2.11) 


after having made use of equation (1.14). When equation (1.15) holds, this 


may be written as 


if tp 
25, =v V Enapy — ce (Siy Gs — Sip G)- (2.12) 


It is immediately evident that if 
=O, al S=conseine, (2.13) 
then 2,,,=0. That is, if the flow is irrotational and isentropic, then 2Q5,=0. 


The converse also follows; for, if the right-hand side of equation (2.11) vanishes, 
we may multiply by U” and sum to obtain 


("is") Sie=0 


in view of equation (1.8). Hence, 
ve US 


Mapy 


ue ANU 


= 0. 
This implies that 
However, it is a consequence of the definition of v" (cf. equation (2.5)) that 

v'U,=0. (2.14) 


Hence, the first of equations (2.13) also holds. Thus, equations (2.13) are the 
necessary and sufficient conditions for 


3. Algebraic Properties of 2,, 
Associated with the tensor 2,,, there is the dual tensor defined by the equation 


Opa Ere), (3.1) 
It is a consequence of equations (2.8) and (2.12) that 


Cl — (ver or +2 su, Bere) (3.2) 
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Hence Q, Nene ee ble a ee (3.3) 
Cc 


Thus, the tensor @,,, will be singular if, and only if, v‘ TS ,=0. If it is singular, 
it can only be of rank zero or two. 
Another invariant associated with the tensor 2,,, is 


2,= 0,0" =2|e* oa, +2 S567). (3.4) 


Since both v“ and TS , are space-like vectors (cf. equations (2.15) and (1.8)), 
the quantity 2, satisfies 

Q, <0. (3.5) 
The equality holds if, and only if, 
v,=0, and S,=0, 
that is, if, and only if, 


It is a consequence of equations (2.11), (1.2), and (1.8) that 


Po an fi 
2,,,U ee Ss 
v de v 
2,0 erro G oP (3.6) 
If av oes to 4 
Quy 8,08 =u, — (eS. S.0)U,, 
where 
w, = E,,910° V7 = gt S 7 
7 UVvOT C2 g ne (3 ) 
and is a vector orthogonal to u“, v‘ and TS _,. It further follows that 
os 26 |(0 2. ay T 
Qyy W =e (2 2) S| Ca (vq, v ea S,,| (3.8) 
and that 
vy —2 (oa de? aly ? 
w, w =e *®? | (y* y,) (= Se?"S,)-(GvS,,) (3.9) 


In case 2,=0, that is, Q,, is singular but not all of its components vanish 
(v“ and T'S , are not both zero), it follows from the above that 


Do == 0 
and . 
QV) 
where 
yt = (u a eee wl 
oe (3.10) 


That is, the vectors v’ and y” are annulled by the matrix of coefficients of the 

tensor 2,,,. The vector y“ reduces to U“ whenever w vanishes, that is, whenever 

digs 
c2 


Ve—A 


We 
So: 
where A is an arbitrary scalar. 
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Since 4 
yy, =14+ 2" _w, w! 


(U% Va)? 


and w,w"=0 as follows from equation (3.9) and the application of SCHWARz’s 
inequality to the space-like vectors involved, the vector y” is always time-like. 


4. Tubes in Space-Time 
If z“(x) is a vector field, the solutions of the ordinary differential equations 


ad x 
aw 


=F) (4.1) 


define a three parameter family of curves passing through a hypersurface 2 
in space-time. We shall assume that the hypersurface may be represented by 
the equation 
at — f(x, x2, 24s) (4.2) 
where 71, 72, 7, and 74 are some permutation of the numbers 1 2 3 4. 
We may then write the solutions of equations (4.1) as 


Ma (En) (1 = 1, tg, 15), (4.3) 
where the €’ are the initial values of the x‘, that is 
ee Oe 


and 


is (&', 0) ae fx (&, 0), x (ee 0), xis (E, 0)] P 
If in the surface 2’ we have a closed curve described by the equations 


Ges Gly Tes CSS 
with 
O(%) =S4e5), 
then the equations 
aM = x (E(t), W] = x* (zt, W) (4.4) 


define a two-dimensional surface in space-time which we shall call a tube. The 
closed curves obtained by setting constant values for W into equation (4.4) are 
images of the closed curve on the initial surface. When in equations (4.4) T is 
constant, the equations define a curve which is a solution of equation (4.1). 


The tangent vector to the curves of parameter Tt is given by 


axl 
Vis ae (4.5) 


where W is kept fixed in the differentiation. We have 
CE ee azH (4.6) 


OW SCWCD 


When the vector field z“(x) is set equal to the four-velocity U“(x), W is 
the proper-time and the curves of parameter W are called the world-lines of 
particles making up the fluid. The four variables &', W may be interpreted as 
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a new set of coordinates in space-time, co-moving coordinates. They are the 
analogues of Lagrange coordinates in classical hydrodynamics. The curves of 
parameter t (W fixed) are then said to “travel with the fluid”. 

Another example of a tube in space-time is obtained by setting z“(x) equal 
to the vorticity vector v(x). Then the curves of parameter W given by equation 
(4.3) are said to be vortex-lines, and the tube is called a vortex tube. 


5. Invariant Integrals 
Consider the integral 
TW) =SV, (x(t, W)] adr, (5.1) 
where V, is given by equation (2.2) and is considered as a function of t and W 
in view of equation (4.4) and A4* is the tangent vector to a curve of parameter 
W given by equation (4.5). The integral [(W) represents the integral of the 
vector V, which is proportional to the velocity vector, around a closed curve on 
the tube defined by the solutions of equation (4.1). 
We now compute 
ryan 
dw i 


t% 


ao 4% ale 
BMA, ede 
This may be written as 


aT 
dw 


Te 


(Q,92° At + 28 Me + Vga) 


= f (Grose + 24 Vi, oa) ae 


where we have made use of equation (4.6). Thus 


“= "|e +o, #) dt. 


Ty 


Since we are integrating around a closed curve, we finally obtain 
ihe 
Sav = | Peot Mar. (5.2) 


The necessary and sufficient condition for /’(W)=0 for tubes defined by the 


vector field 2“ by means of arbitrary initial surfaces and arbitrary curves in 
this surface is that 


Qa 8 — O . (5.3) 
Similarly, the necessary and sufficient condition that d/'/dW =0 for tubes defined 
with the same amount of generality is that 
as B ead Zp, a = 6) 
where 


24 = 245%. (5.4) 
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6. The Circulation 


In this section we set 2*= U%, the four-velocity vector and W=s, the proper 
time. We shall then name the integral in equation (5.1) the circulation and 
denote it by C(s). Equation (5.2) then becomes 


dC * 5 Op a 
fee 35 Sah at (6.1) 
in view of the first of equations (3.6). 


Thus, a necessary and sufficient condition for the circulation to be inde- 
pendent of proper-time is that 


~~ 9” _ (,-~ 9P 
(e ae (e f S;4) (6.2) 
It is clear that for isentropic flows (S=constant), dl /ds=0. 

When equation (1.15) holds, equation (6.2) may be written as 


(LS, .);8 — (LS;¢),.=0. (6.3) 

In virtue of equation (1.9) this may be written as 
P;a0;8 — P;60;a =O. (6.4) 
Hence, in this case the necessary and sufficient condition that dC/ds=0 is 


that at every point of space-time TOP (6.5) 
That is, the motion is such that a knowledge of # (or g) at a point determines 
@ (or p). 

7. Bernoulli’s Theorem (Time-dependent Flows) 


The definition of the circulation, namely, 


CRS) = (reG UA dx 
involves an integrand which, aside from the exponential factor, is the tangential 
component of the four-velocity field U* taken at a point of a curve which travels 
with the fluid. C(s) is, thus, similar to the classical circulation. To show that 
it is a generalization of this classical quantity, we shall show that as in classical 
hydrodynamics, the condition Gisk==0 (7.4) 
and is related to Bernoulli’s theorem for time-dependent irrotational and isen- 
tropic flows. 
It follows from equation (5.3) that the condition (7.1) is equivalent to 


2.,,=0, 


that is, to the requirement that the flow be isentropic and irrotational (cf. Sec- 
tion 2). Under these conditions (cf. [3]) there exists a scalar function R(x) 
such that 4 oR 

Vi=e*U, = (7.2) 


Oxt © 
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Define the quantity 


(dh MTT gore a0 ( fan bs») U#U’. 
eke bur Canaub 844 ese 844 
This is the length of the vector 0” which is given by 
(Ola —_ U" {a4 Of U, 
844 


and which is the projection of U“ onto the three-space orthogonal to the curves 
of parameter x4. 
In view of equation (1.2) we have 


I 


eee || 2 
U,=Veas 1+. (7.3) 
Equations (7.2) with w= 4 imply that 


peel Vi Read kn 
e ?’U,=e "Venu //1+t= 


ox © 


When equation (1.15) holds, this becomes 
) a 2 OR 
(1+ 5) Vea 1+ $=. 


2V 
Sage uae sao (7.4) 


If we now write 


and expand the left-hand side of the above equation neglecting powers of 1/c 
higher than the second, we obtain 


(i HV +4 ?) = 2+ ee 1). (7.5) 


ox* 


This is the classical form of the Bernoulli equation for a non-steady irrotational, 
isentropic flow in the gravitational field with potential V. The quantity c?R 
is to be interpreted in terms of the velocity potential for the flow. 

Equation (7.4) is the usual relation between the classical Newtonian potential 
and the metric tensor of space-time which holds in the case of weak fields. 

In case we use equation (1.13), equation (7.5) becomes 


(+0 +54) =le+ilby, )\(Z% —1). (7.6) 


On 


8. Bernoulli’s Theorem (Steady Flows) 


In classical theory Bernoulli’s theorem is usually formulated for stationary 
flows. LICHNEROWICZ [2] has indicated a method of describing steady flows of 
fluid in their own gravitational fields. It consists of assuming that the space- 
time and the velocity field are invariant under a time-like one-parameter group 
of motions. That is, there exists a vector field é* which is such that Killing’s 
equations 

ooh pa Sty =O (8.1) 
Ue UO, 
Sas Pine = 0s 


are satisfied and such that 


(8.2) 
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Equations (8.2) state that the velocity field U’ is invariant under the group 
of motions as are the entropy and pressure. 
Define the quantity 

Hese 2U eae (8.3) 

and compute 
fey ees (Cmca, Us 
We have 
Hg U? = e-*(— gp, , UPU,& + U,, UPS + U, Ey UP) 


= 60 | Ba + UFUP Eg + by: 0)| =0 
as a consequence of equations (1.7), (1.8), (1.14), (8.1) and (8.2). Thus, along 
a world-line of the fluid we have 

H = constant. 


If the coordinate system is chosen so that 


i= of 
this equation becomes 
én*-U,=‘constant: 


Using the results of the preceding section, we may write this as 
a+ oe + V =constant (8.4) 


along a world-line of the fluid. It may very from world-line to world-line. Equa- 
tion (8.4) is the usual form of the Bernoulli equation. 
It is instructive to compute the derivative of H in an arbitrary direction. 
Thus, 
Hg =e *(— 9p Uo + Uz; + ULF 6] 
= ole Qa gS — Pa 5" Us + Ups aS — UNE g + U* (bai +85; a) 


as follows from the definition of Q, , (cf. equation (3.2)). In view of equations 
(8.2) we have 

A, 5 = 2,55. 
Hence, if 2,,=0, 

H = constant 


thoughout space-time. That is, if the motion is isentropic and irrotational, the 
Bernoulli constant is independent of the world-line of the fluid particle as is 
the case in classical theory. 


9. Vortex Tubes 


In this section we shall discuss the integral (5.1) when the tubes in question 
are vortex tubes, that is, when we set 2“(x)=v"(x), and the motion is rotational 
(v‘=+-0). In this case, we cannot have /’(W)=0. However, it follows from 
equations (5.2) and (3.6) that 


ih rs in OP » a 
Mee aie 35 5.0% Cintas (9.1) 
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Equation (9.1) may be written as 


dT 
a ee 4 
Sg =f Ved dt (9.2) 
with ag . 
aoe (9.3) 


The necessary and sufficient condition for dl'/dW =0 for tubes defined by arbitrary 
curves in these surfaces is, then, 


(Vz); p — (V;);2 = 0. 
This may be written as 
hQugth pV_—haVp=0.- 
Thus, 
hQ, p= aVe— hp Ve, 


and is of rank 2. Therefore, h=0O or 2,=0. However, £22,=e°h. Hence, the 
necessary and sufficient condition for d//dW=0 for vortex tubes defined as 
above is that 2,=0, and Q,, be of rank 2. 


In this case we also have 


dT 
rier 


for world tubes defined by the vector field y” given by equation (3.10). 


10. The Change in Circulation Across Shocks 
All of the preceding discussion was based on the assumption that the deri- 
vatives of the flow variables existed everywhere. We now consider the possi- 


bility of singular hyper-surfaces in space-time, representing shocks across which 
equations (1.6) obtain. These equations may be written as 


™ = 0, UNA,=@ UMA, (10.4 
and 
AE 
= (u, Ut —p_ UH) == (2, —P), (10.2) 
where 
a rebag Ute ei (10.3) 


In case equation (1.15) holds, equation (10.2) may be written as 


m ; Ae 
* (Ve Vy =< (b. — f), (10.4) 
where 4“ is normalized so that 
MA, =—A (10.5) 
and for a shock wave m=+0. 
It then follows from these equations that the components of the vector 
field Vy, tangent to the hyper-surface are continuous across it. As has been shown 
earlier, the components normal to the hyper-surface are discontinuous. We may 
determine the discontinuity in this component in terms of the discontinuity in 
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pressure (or density) across the shock as follows: multiply equations (10.2) by 
U,,,,U_, and A, to obtain 


=p 


(#, — #_UEU,,) = ——(p, — B.), (10.6) 
(#, UU, —w) =< (6. —B-), (10.7) 
Pilg Jotaes eg ae (10.8 


respectively. If we now eliminate U“U,,, from the first two of these equations 
we have the equation 


1 A : 
M-w@=3 6-2) (+4). (10.9) 
From equations (10.7) and (10.1) we have 
Ue PsP en\? Q+ | 20= 
ee (aa ai eA) Se WEA) - (10.10) 
Q- Q+ 


Equations (10.9) and (10.10) are the relativistic Rankine-Hugoniot equations 
previously derived [4]. When yw is known as a function of # and 9, equation (10.9) 
may be used to determine #, in terms of #_, g, and 9. Equations (10.10) may 
then be used to express 9, and hence #, as a function of V“4,, p and e_. They 
may also be used to express V/A, as a function of the latter quantities. 


We shall write 


Nees een ee a yet 10.14 
(oot Teepe Mody) 
Then equations (10.9) and (10.10) may be written as 


M—1=(1— =) (Vea) =o8(y— 1) (+4), 


+ “ue 
(10.12) 
pe sate fey 
y lo 
pane 2 
cep” 


When yw is known as a function of # and 9, the first of these equations may be 
used to determine y in terms of 7. We may then use the remaining equations 
to determine all quantities as functions of V"4,. Since the entropy will not 
be conserved across shocks, care must be taken to choose solutions of these 
equations such that 

Sp Sz 
across the singular hyper-surface. 

When a shock is such that 4, is a function of the points of space-time, we 
say that it is of varying strength. This variation may be due to the curvature 
of the shock at each time x4= constant or to the dependence of 4, on x* alone. 
The world-lines of particles which cut the singular hyper-surface have discon- 
tinuous tangents at the hyper-surface. Along such a world-line the entropy 

22* 
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suffers a discontinuity. The size of this discontinuity in entropy varies from 
world-line to world-line. As a result, it is not true that the pressure and density 
at all points of space-time are related as in equation (6.5). Thus, when shocks 
are present, we cannot have dC/ds=0 for arbitrary tubes. 

This result may be seen in another manner. Suppose we have a shock of 
varying strength (A“ not constant) and suppose it is such that V", p_, and 9_ 
are constant. Then, we would have 

Qiy =F Vee rr ane = 0. 
However if 
Qi =0 , 
then it follows from results of Section 2 that the flow behind the shock would 
be both isentropic and irrotational. However, we have seen that it is not isen- 


tropic. Thus, 
Shek = AY 


This means that we can construct tubes of world-lines which cross shocks and 
for which C(S)=0 on one portion of the tube, but C(S)=-0 on other portions. 

The detailed behavior of the circulation across shocks depends on the explicit 
expression obtained when the system of equations (10.12) are solved as described 
above. These expressions depend in turn on the nature of w as a function of # 
and @ or equivalently on the nature of ¢ as a function of # and g. Once this 
function is prescribed, we may study the geometry of the world-lines in terms 
of the geometry of the shock surface by introducing a coordinate system con- 
sisting of the world-lines of the fluid particles as curves of parameter x* and 
the singular hyper-surface as the hyper-surface x4=constant. The methods 
used in studying stationary and pseudo-stationary flows [5] may then be applied 
in a straightforward manner. 


This work was supported in part by the National Science Foundation. 
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On the Convergence 
of the Rayleigh Quotient Iteration for the Computation 
of the Characteristic Roots and Vectors. II 


(Generalized Rayleigh Quotient and Characteristic Roots 
with Linear Elementary Divisors ) 


A. M. OSTROWSKI 


In this part * of our paper we give a generalization of the Rayleigh quotient 
iteration method which can be applied in most cases to non-symmetric matrices 
(in the complex case, to non-Hermitian matrices). The enumeration of the 
sections and formulae continues that of PartsI and Ik. }, 34; 


25. The reason why the Rayleigh quotient method works out so well in the 
symmetric case is to be found in its connection with certain variational problems, 
since in the symmetric case the eigenvalues are the values of the extrema of the 
usual Rayleigh quotient, corresponding to the eigenvectors. 

In the case of non-symmetric matrices, however, the quadratic form in 
question has to be replaced by a bilinear form. But in this case no extremal 
property analogous to the former one is known. However, asa matter of fact, 
the reason for the good convergence of the Rayleigh quotient iteration lies not 
in the extremal property itself but in the stationary property of the Rayleigh 
quotient in the neighborhood of an eigenvector. Under these circumstances it 
can be asked whether the stationary property still holds in the case of a generalized 
Rayleigh quotient, and this turns out indeed to be the fact under fairly general 
conditions. In the corresponding cases the generalized Rayleigh quotient iteration 
remains convergent also for non-symmetric matrices; the convergence turns out 
to be cubic, and even the main result (46) of the symmetric case remains generally 
true in the non-symmetric case. 


26. For a general square matrix, A, of order , we shall call the expression 


(72) R(En) = 


the generalized Rayleigh quotient (to be denoted by G.R.Q.) corresponding to the 
vectors &€ and 7. 


* This work was performed in part at Numerical Analysis Research, University of 
California at Los Angeles, and in part under a National Bureau of Standards contract 
with The American University, Washington, D.C., with the sponsorship of the Office 
of Naval Research. I am indebted to Mr. G. CuLLER and Dr. E. V. HAyNsworTH 
for discussions. 
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We assume through the end of Section 28 that A is real. Corresponding to 
a real eigenvalue o of A, let &) and % be (real) right and left eigenvectors of A: 


(73) A&=o8, oA =o%- 
Assume further that 

(74) né'+0. 

Then we verify at once that 

(75) R (Eo, No) = 9: 


We ask now whether R(&, 7) has at the point (£ , 7) the stationary property 
with respect to (&, 7), that is to say, whether the partial derivatives of R(&, 7) 
with respect to the components of € and 7 vanish at the point (&), 7). We shall 
show that this is indeed the case if o is a single or multiple eigenvalue of A, 
provided that all corresponding elementary divisors are linear. On the other hand, 
we shall show by an example that if a non-linear elementary divisor corresponds 
to o, then the stationary property need not be satisfied. 


27, Assume that o is an eigenvalue of A of multiplicity m with only linear 
elementary divisors. Then we can introduce special coordinates in the &-space 
and in the 7-space in such a way that the first m coordinate vectors in the -space 
become a complete set of m linearly independent eigenvectors of A from the 
right and that the same property holds also for the 7-space. Then, if the com- 
ponents of € and 7 are correspondingly denoted by (%,,..., x,,) and (1, .--, Y,), 
we have 

nAF =a(% y+ ra: +n Ven) + B(Xmii> very Xny Vmtarcees Vu) 


where B(%y11, +++) %n3 YVm+1» +++» Y,) is a bilinear form in the variables x, and y, 
written out as arguments. 


Now let € and 7 tend to the eigenvectors &, 7, where 
Sq (Gy, 1p Ee ORte 0) Fo =O OLS TOROS 
Grae Ay 5 aus yg oh Coca > nun opal) an ment On eae nh ener a 
and 19 &) = 40, +++: +4 ,b_=- 0. Suppose further that we have for e+>0, %,=O(e), 
Bu.=O(e) (u=m+1,...,n). Then we have 
nA&' =a(a,b,+---+4,,6,,) +0 (e), 
n&' =a,b,+---+a4,,b,, + O(e), 
Rin) —0 =) 


fr 


nS 
WE >No b= 0; 


and since 
we have finally 
(76) R(E,n) =o + O(c). 


This result is the indication of the stationary property in this case. To give 
the formal proof, observe that we have 


Aes ) —~o= Blom+a, vee, Oy; Batis ses» Br) = (Cpa, Bara i Op Bn) : 
(4,5, 0+ + Am Dm) + (m+4 Bmtits ten Bn) j 
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but here the numerator on the right is of second order with respect to the co- 
ordinates, which vanish in (&), No), While the denominator tends to m &j. This 
proves our first assertion. 


28. To prove the second one, consider the matrix 


Dae ( o 4! 
ONG; 
which has o as a double eigenvalue and the elementary divisor (A—o)?. The 
eigenvectors are from the right (1, 0) and from the left (0, 1). Here, if & =(x,, %9), 


Nis) thea. iO) is 
0 (4%, Vy +%2 V2) +% Vo ; 
44 Vit %_ Vo 


For &=(1, a), 7 =(8, 1) with «, B tending to 0, this becomes 


1 
aR rer try i 


and there is no trace of the stationary property at the point &)= (1,0), 7)>=(0, 1). 
29. The generalized Rayleigh quotient iteration corresponding to G.R.Q. 
defined by (72) is given by the formulae 
(77) Ay = R(x, Ne)» 
(A — 4,1) fy. =&, 
Mold — Al) =m, 


which are to be applied consecutively for k=0,1,..., starting with a couple 
&),%- These starting vectors are chosen arbitrarily, subject to the condition 


No £0 0%. 

Let o be an eigenvalue of A to which correspond only linear elementary divisors. 
Then we are going to prove: 

Theorem 1. There exist a positive )=6(&, 49, A) depending on €,, no, and A, 


and a certain algebraic function of the components of the vectors € and n, p,(&, 7), 
depending only on A and the choice of o, such that 1f we have 


(78) 


(79) [Ay — o| S 6, 
(80) P1(So, %0) + 9, 


then the sequence A, tends to o. 


* Another approach to the generalized Rayleigh quotient iteration has been given 
by M. R. HesteEnes, Inversion of matrices by biorthogonalization and related results, 
J. Soc. Indust. Appl. Math. 6, 51—90 (1958), who derived it by applying Newton’s 
method of differential corrections simultaneously to eigenvalues and to the correspond- 
ing couples of eigenvectors. From this point of view the cubic convergence of our 
method is particularly interesting, since Newton’s method usually gives only quadratic 
convergence [see, e.g., A. OSTROWSKI, Konvergenzdiskussion und Fehlerabschatzung 
fiir die Newtonsche Methode bei Gleichungssystemen. Comm. Math. Helv. 9, 79—103 
(1936 37) and A. OstrowskI, Simultaneous systems of equations. National Bureau 
of Standards, Appl. Math. Ser. 29, 29—34 (1953)]. 
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Theorem 2. If the sequence A, tends to o, without any of the i, becoming =o, 
and &,, 1 satisfy the condition (80), then we have 


(81) Ay+1 — 0 = 0[(A, — 9)?] (k > 00); 

moreover, for any fixed integer Pp, 

(82) Ania — Oo = OLA, — 0)" (Ay ey, 40) a (Age 8) (k > 00). 
Theorem 3. Suppose that under the hypotheses of Theorem 2, one of the two 


following conditions 1s satisfied : 


A. Denote by o, all eigenvalues of A different from o and different from each 
other. Consider among these o, those whose distance from o 1s minimal and denote 
them by 01,0,-..,0,. Let L be the greatest exponent of an elementary divisor 
corresponding t0 0,,09,...,0,, and assume that to only one among 01, 02, .--, 0; 
really corresponds an elementary divisor with exponent L*. Further, assume that 
a certain algebraic function of the components of £5 and no, p3(§o, No), depending 
only on A, &, Ho and oa, does not vanish. 


B. All eigenvalues of A are real, and A has only linear elementary divisors. 


Then we have 
(33) MHF 5g 40 (ko). 


30. Assume that o is a characteristic root of exact multiplicity m=1 and 
denote the elementary divisors of the matrix A corresponding to the eigenvalues 
=-o, in a convenient order, by 


(84) (A — 0,)",..., (A — o)™, 


l 
(85) mt>in=n, Maxn,=N. 
v=1 ‘e 
If A is transformed by a matrix S to the Jordan canonical form, we obtain 
for S14 S the representation 


\ 1 
(86) oO ay + >: T,, , 
= 
(87) T= 9, Ls ae es 


where generally J, denotes the unity matrix of order k and U, the so-called 
auxiliary unity matrix of order k consisting of units in the first superdiagonal 
and zeros elsewhere. The symbols +, > denote the direct sum of matrices 
aligned along the main diagonal. In this transformation the vectors MH, are to 
be replaced by »,S™, and the vectors & by S&,. The inner products 7,& are 
then, of course, invariant, and the same holds for the generalized Rayleigh 


quotient defined by (72), and for the equations (78), and of course also for the 
values of A, given by (77). 


We can therefore assume from the beginning that A has the form (86). 


i Of course this does not exclude the possibility that this elementary divisor is 
multiple. 
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31. Let the decomposition of the vectors &, and , Corresponding to the 
decomposition (86) of the matrix A be 


(88) APD, naa ta 
and put 
(89) b= p= Pi (So. Mo) = Ne ue 
B=>>T,, 
(00 cai 
Eee cay — ee, 
y=1 v— 


Then we have from (78) 


1 1 
(91) ElEY 7 a EM), nee = =n? (k =0,1,.. .) , 


02) (BA, Nee Hs, nf BAD =n — (k=0,1,...) 
( 7 1 


From (91), putting 


k 
(93) Me Li ic 4) (R= O; 15 <2); 
we have 
() ee (0) (2) eat ITE} 
(94) gy = Mees ei, Gi = Vn m1: 


(95) mA = ame 42, 7) — pl) Bee, 
k—- 
(96) “neh = vale eNO NOE nf elt)’ 
—1 
and from (77), (72) 
Z(*) 6 N(&) 
(97) a le oes - = Mz-1D,, 
Set EN 
Mente 
Z*) _o N(@*) 
= Pe 2p, + M2_, N®)* 
Now put 
(99) Neate, 20 ter, 


Then from (90) and (92) we have, since U,, commutes with E,,, 


Demet alers. (Ss Nase lag tel y,) = Thy 
k-1 k— 
ge — Cr oe Ae Hoe ie EP)” ) ne) = ne, i, (le, A, I ,,) 2 
o—0 x=0 
k—-1 f 
(100) Te yee be das) 2 Sd, ye 
x=0 
k-1 ’ 
(101) Ze Le eS WL No en: 
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From (90) and the definitions of Z and N™ we have now by (99) 
l l 
(102) 2 =S 2, NY =P NM, Zo NO => (Zo NM). 
y=1 y=1 


Further, from (97) and the definition of M, in (93) 


Inti —F _ Dati 
WD (A,—o)? ~~ Dy 


33. In what follows we shall have to use the notion of the bound of a matrix. 
The bound A(C) of an m xn matrix C is defined by 


(104) A(C); SanaxhCF| sale ead). 


From this definition, we see at once that if the matrix C is majorized by a 
matrix C, with non-negative elements, C<C,, that is to say, if the modulus of 
each element of C does not exceed the corresponding element of C,, we have 


(105) A(C) S A(G) (C<C,). 
From (104) follows the triangular inequality, 

(106) A(C, + Cy) S A(Gy) + A(C2) 

and the homogeneity for any scalar, c, 

(107) A(éC)=Tel Ale), 

as well as 

(108) A(C3C.) = Ale, JAG.) 


It is well known and easily proved that A(C)? is the greatest eigenvalue of 
the matrix CC*. From this it is seen at once that A(£) =1 and it can be easily 
deduced for the auxiliary matrix U, by forming UU’, that A(U) =1*. We obtain 
now, for any scalar, w, 


(109) A(jwI + U)SA(wl) + A(U) S|w| +1. 


Further, by using (105), (106), (107), for any scalar w, with |w|=1, since 
Ue =O; 


1 1 U U2 Ur-1 d 
Alereu)— Torts tot 7 Se) a 
(110) A( oa )sn (|w| 4). 
34. Now put 
(111) Min |, — o| oy) 


The number d is obviously a homogeneous function of the elements of A in the 
sense that if all elements of A are multiplied by the scalar c, then d is multiplied 
by |c|, while o and o, are multiplied by c, and the same is true of the A, deduced 


* Cf. Ostrowsk1, A. M.: On the bounds of a one-parametric family of matrices. 
Journal f. d. r. u. a. Mathematik 200, 190—199 (1958) 
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from (77) and (72). The three assertions stated in Section 29 remain invariant 
if all elements of A are multiplied by c=-0, though of course the values of the 
constants 6 and g occurring in these assertions are changed. We can therefore 
without loss of generality assume that 


(112) |o,—o|=2 (o, + 0). 
Further, put 
l l 
(113) Ps = fo(Eo, No) yi | 
v=1 y=1 
(114) p= (Eo. = PEL, 


Now we make the assumption that the vectors &, 7), though arbitrarily chosen, 
satisfy the inequality 


(115) pele 0. 


If one of the vectors &, 7% is already an eigenvector corresponding to o, 
then we have #,=0, f= oo. Therefore the size of f is to a certain degree an 
indication of the nearness of one of the vectors &), 7) to a corresponding eigen- 
vector. On the other hand, if one of the vectors £), 7) had only vanishing com- 
ponents corresponding to the eigenvalue o, we should have #=0, and this is 
excluded by (115). 


35. Now put 
(116) Max|o,—o|=4—1. 
From (100) and (101) we have 
(117 29) — NP =O Tay — 0nd) TT (Tay — Pad) 282 
Here, by (109) and (87), we have obviously 
(118) A(T, —o1,,) S|o,-—o0| +18 4,. 


We now make the assumption that for a fixed value k considered in this 
discussion we have 


(119) |4,—o| 1 = 04,...0h —1). 
Then from (112) follows 
(120) |o, —A,| 21 (octe Onl ok aoe i aD eae da 


Using this we obtain from (110) and (85) 


1 1 oe AGS Foe 
(121) MN reg | = eyiaar) Na (c=0) 450, PN =A). 
Using (108), we now obtain 
k-1 
(122) Al (Lay — 0 Fy) LT Tay = Be Toy) | S 4, N* 
x=0 
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From the definition of A and (117) we have 


(123) |Z? — o Ny] S dyn" | ney | Ee 
and, further, trom (102) and (113), 
l 
(124) 12) — GN] < dy Yn * |) [£4] < dy ON 
y=1 
Moreover, again using (121) and (110) we have 
k 
(125) Al TT (Tay — Ae Zne)3] S02" 
“x=1 
and therefore by (100), (102) and (113) 
(126) [NS] <0f* | PEO] = ny" NO, 
l 
(127) JN | <i n2* [NL] < pO) NPA, 
vy=1 
36. We have now from (98) by (124) and (127) 
d iP N2k 
sos Pal <= OPN MP 
Define 6 by the formula 
Min (3g. 
(129) 6 = 0(€ 9%) = — aa A 


and assume, instead of (119), that the sharper relation 


(130) lo—A,|<S6 (= 0, .: SRE=4) 
is verified. | 
Then by (93) MAS of, 
and we obtain from (128) by (414) 
\D,| ze a, po N2k = a,N™* i < a,N** 
2|p|— pO (SN)2® 2b —dN2. 6-(dN)2?*-® — 2p—dN? 


It follows, since d,=1 and therefore by (129) 6N?<4, that 


|D,| <n. 
We see that by (97) 
(131) Sais SBN" IMP 


and therefore 


|4,—o] <2 (ON) 


~/5 


= on 4 SN") (ayy 1 gh. 


But here in the right-hand expression by (129) the second and third factors 
are both <1, and therefore we obtain 


(132) A= oO S08 
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From (130) it follows, therefore, that the same inequality is also satisfied for A,. 
Therefore (132) remains true if we assume instead of (130) only 


(133) |o—A| <6. 
37. In order to improve (132), we note that by (132) 
(134) |M,1|? 5 IL Oye Sralin 
and therefore from (131) rs 
Was (S SOONG aaa 
here, by (129), the factors in parentheses on the right are <1. Thus we obtain 
(135) Are Ole eet ee), 
Putting k +1 instead of k in (131) and dividing both sides by the product 
(A, —o)? (A,_1 — 0)... (A,_; — 9)?, 
for a fixed ¢<k and by the definition of M, we obtain obviously 


Anti — 9 


t 
II (Ag2= a)? 
t=0 


d 
(136) Ss pe | Mal? 


The expression on the right, by the definition of M, and by (135), is 


- 4, N2k+2 | AS of sonett guerre a (6.N2)* HL g2h*—(4t+9)k4 142 (44+2)* 


and by (129) this is 
< 2B (t+ 9) +2 (42/41 


and tends to zero with k— oo. 

The convergence condition (133) with the expression (129) for 6 has been 
derived under the assumption that d as given by (111) has the value 1. If we 
now drop this assumption, then the general case d=1 is reduced to the case d =1 
by replacing A,, 0,0, and d, by a + cr, BEES The condition (133) 
then becomes 


fear ray pd 
Se 
< Min ieee 


or 


and we see that in the general case the number 6 used in (79) is defined by 


ik ec abe tg Ms 
(137) jo = 4 Mi ee 5 (Ey, 19,0, A) 


The initial vectors &),7 enter into the bound 6 given by (137) only in the 
expression p defined by (113) and (114). 
Our Theorem 1 of Section 29 is now completely proved. 
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It may be observed finally that if the number N defined by (85) is not known, 
it can be replaced by », the order of A. On the other hand, if A has only linear 
elementary divisors, we have N =1. 

The iteration defined by (77) and (78) starts with a couple of vectors &), No» 
and condition (79) implies that &) and 7 are a more or less good approximation 
to certain eigenvectors of A corresponding to o. Indeed, in this case the formula- 
tion of the stationary property discussed in Section 27 shows that Ag will be 
then a good approximation to o and condition (79) is likely to be satisfied. 

If, however, we wish to start with an approximation to ao, we denote it by 
A_, and from (78) applied to k = —1, derive the values of & and No» taking for 
&, and 7_, arbitrary vectors. Then, as we have already discussed in a previous 
publication, the vectors é) and 7, can be expected to be good approximations 
to certain eigenvectors of A and our discussion will apply. 


38. We are now going to prove Theorem 2. In fact, the assertions (81) and 
(82) of this theorem have already been proved by formula (136), but only under 
the conditions of Theorem 1 and not under those of Theorem 2. It is therefore 
sufficient, in order to prove Theorem 2, to show that if A,—>o and (80) is satisfied, 
then for a certain k we have 


(138) [Ap o| =O = Or ne 4) 
and 
(139) Pils Me) + O- 
But we have from ee and (94) 
vs No 
(140) = (Ess) = al £0 = Paka tel 


and since by (80) the numerator does not vanish, (139) is certainly satisfied. 
In order to discuss 6(&,,7,), observe that by (113) and (126) we have 


(141) ps = pa(&,. 2) = 2 IN| sa DN] =n 
and by (114) and (440) 
(142) bp” = b(&.m) = Pools woke [PL 


pe» = oP VERS po 

Here, since 4,0 (k-> 00), we have n?* M?_,->0; therefore 
(k) a it Git 5 i Pp) d | 
60 = 8 Esme) = aye Min (>, Peg 


from a certain k on becomes equal to d/(2N2), and therefore for a certain k = Ro 
we surely have 


(143) lo — A, | < 6) 


But then if we put 2,=A,,,,, the assumptions under which Theorem 2 has been 
proved are satisfied for A,, and the assertions of Theorem 2, being true for A,, 


then follow immediately also for the sequence A,. Theorem 2 has now been 
proved *. 


* An analogous argument has to be added in Section 17 of Part I of this series, 
where the corresponding part of the proof was not developed in all necessary detail. 
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39. In proving Theorem 3 of Section 29 under the hypothesis that A,—o, 
we can again assume, without loss of generality, that (79) is true. 


Before discussing the limit in (83), we prove the following 


Lemma. I7 I and U are the unity matrix and the auxiliary unity matrix of 
order n, and if a, 1s an infinite sequence tending to a=~0 as k-> 0, then 


k 
1 1 Oe n—1 
(144) wo TT Ta > Cane U (k + oo). 
Proof. For sufficiently small positive x, put 
k 1 lee) (®) 
es fe 
(145) Hq = > se, 
x=1 yu=0 
where Svea (a,.--,@,) iS a symmetric polynomial in a,,...,a,, homo- 
geneous of degree w. Then obviously, since J and U commute and U”=0 
k n—1 
1 ad (k) 
eer Po 
“x=1 u=0 


and (144) follows at once if we prove that as kono, 


Sh L 
(146) ae (u=0,1,...,%—1). 
Put 
On b> (oc ed AO at) 
a 


Then (146) is equivalent to 


SOG tts Op) 1 


(147) Re < be! 


(k + 00). 


Therefore from the beginning we can assume a=1, a,—>1 and of course u >0. 


oot saat +: +a —(i+eqk - Y=0Ay2%)s 


where ¢, ,—0 (k—> co). Then we have 


: (R) yo 
Shae 
IT 1, Fie Os yan, 
x=1 v=1 
and therefore 
(148) exp|& >) (it &,4) 2] =1+ 4? 2, 
y=1 y=1 
where 
(149) AOS Vl aA ale) 2 
Putting here x = y/k we obtain further 
: y” Se 
(150) exp{(ttadat D+ 6,0) yporp= t+), ey 


As k-> oo, the expression on the left-hand side tends to e*, and the coefficients 
of its development in powers of x tend to the correspondent coefficients of e”. 
We obtain A®/k’—1/y!, and (147) follows from (149). 
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41. From (117) and (87) we have 


k—1 


Y) — @ NIP = 98), ((ay — 2) Let Ung) ! ((0,::Ag) Lapis) ct Barer 
a 
(151) 22 —o NY = 98), ((0, — 6) Lup Une) TL aoe Ha Sa. 
where 
(152) TT, = if (In, +5 Un) 
But here we have 
Oeuaa = Cine a a> ae =a, 


and therefore our lemma can be applied to the product //,. Observing that the 
number of factors in //, is 2k, from our lemma we have 


IT; (CEA cat. fy —1 
(Rye s Gea U (k + 00). 
Since U*”—0 as k— oo, we have from (151) 
Zo N® ise —2)""—-1(6,—0)-™71 EAs ey, 
(153) ee an (a, + Ae) > ( ) (n (6 1)! (7, Bb : E; 0) ‘ 
x=0 a s 


For the product on the left-hand side of (153) we have 


k—-1 L k—-1 ae ‘ 
IT © —4) =(,— 0) J] (1+) =o} 2, 
x=0 “x=0 % 
where 
rs 
Te ee ps4) 
Therefore from (153) we have 
(Rk) (k) —_2)"— 
(155) Lis a (o, om < (—a)teat . (n®, ica &,) (ps oo) : 


[m= 11 


42. We assume now that the conditions under A of Theorem 3 in Section 29 
are satisfied. Then, if the o, in (84) are suitably ordered, we have 


(156) Min|o,—o| =|o,—o|=---=|o,—o|<|o,,,-o|S---<|o, 
(157) Leng =n, > nye SN, 
(158) 01 = 0,=-+- =9,. 
Then from (154) we have 
R=P=--=P, 
and it follows from (155) that 
ty apts, Ya YUP plosegal ppralesey 
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If now we assume that 


(160) b= Dnl, UE e+ 0, 
and this is an algebraic coKaitien imposed upon &, and 7), then we have 
Z) go NYW wy, poh ess ae i) 
(161) hh ls = 
A mea Ceenia a uae 


On the other hand, it follows from (127) and (134), since 6<1, that the 
denominator of D, in (98) is ~nfP é)’ +0, in virtue of (115). Therefore from (98) 
and (161) we obtain 


pl-l 


Y. 
(162) D, es i) 507 te (o,—0)2* (k a oo) F 
Dr4i 1 
(163) = > Gas (k — 00) 


and from (103) we have finally 


ue a herit e eoe 
ea) (4, —a)®* (o,—0)? * 
43. \WWe assume now that the conditions under B of Section 29 are satisfied, 
that is, that o and allo, are real, and N =1. Then, since each n,=1, (117) becomes 
k-1 


(165) Zy) —a Ny = (ns), €8 ) (0, — 0) IT (0, — 4.) 


oO) 


For the product on the right, using the decomposition 


1 eee Fa (1 2— | 
6,—A, 6,—6 G5—A,]’ 


we can write 


k=l 
Il (Gp 1,)? = (6, —0)7** Q, 4, 
(166) ses naa aati 
Q.0= TT (1 i] 


and therefore 
(167) LP ONPG yes eee 0) = Ose: 
From (166) we have 


o—A, \* o—A, { 0,—C 
Q,,241— Ov = O,0|(1 = =) =H = Q, 4+ ( j 1), 


oO — A; Ap—Oy 0,—Ap 
v,k Oy 
(168) [Q,,n41— Ona] S| ee (SEE + 1)| [o — A 


44. Now, if the conditions (79) and (80) of Section 29 are satisfied, we have 
(135), and therefore the series >’ |o—A,| is convergent; further, o,—/, tends to 
ca 
a 


o,—o=0 and the series Ds = 
x es bed 
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is also convergent. We see that the product 
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Q,,, tends to a limit Q, as koe and therefore the modulus of the coefficient 
of lo —j,,| in (168) is not greater than a constant c independent of k andy. From 


(168) we have 


Ope Q 


and it follows that 


|Q, — 2, 2| = PAC v,4t1 Oa 


= 60 an] (1 


from some k on. We can therefore write (167) in the form 


sed 


> ee 


(169) Z{ — o Nf = (mf, &’) (0, — 0)-** 0, + O[(, — 0) -** |o — Al], 


(170) 0,= je —ae 


Denoting the factor (7{), £2’) (6,0) Q, by ¢,, we obtain 


Wen eS 4 of eae! 
a) an aS parent ean 
Put 
1 — 
(172) Max (are ts 
Then by (102) 
I 
(173) ZN = ON ae a | O(|o—A,| 7"). 
y=1 ? 
Denote the different values of the expression ar in decreasing order by 


Then (173) becomes Vanna ri Eo Ye 


(174) ZA GN ie Sepak te nse Spee +s, yf +0(y*|o— A,|)- 


Here some of the s, might vanish. If we retain only those y, to which seats 
non-vanishing s, mad denote them in decreasing order by g, > g, >-:-> g, and 


likewise the corresponding s, by 4, f,...,t,, we obtain finally 
(175) Z® — oN =t,¢8+---+4,e8+ O(y*|o— A), 
unless all of the s, vanish, in which case we would get 
(176) Z" —_ oN =0(7*|o — |). 


45. If we assume that (176) is impossible, then (83) follows easily. Indeed, 
from (127) and (134) (observe that (112) can be assumed to be true) 


[Mpa NM) spe) 62 (N20 )"— 0; 
and we see that the denominator in the expression (98) for D, tends to 2f,--0 
in virtue of (115). Therefore it follows from (98) and (175) that 
Diss 


Lok 
D,~ ees Rate, 
k &1> De Ra 


The Rayleigh Quotient. III 339 


and by (103) we have 


Anti~F _ Dati 
(,—o® Dy 8! 


where g,;=g is one of the numbers ie 
j (o,—a)? 
It remains now to prove that (176) is impossible. Indeed, from (176) it would 
follow by (98), since the denominator tends to 24,, that 


D, =O (y*|o — A,|) 
and by (97) 
A, —9 =O(Mp_sy*|o —4,|), 


; 1 
that is to say, that ae 


would be bounded. But by (134) 
ad 


1 1 1 


ea aay 
would be bounded too, while, since 0<6<1, the denominator here tends to zero. 


Therefore (176) is impossible, our assertion (83) is now proved in the Case B, 
and the proof of Theorem 3 is now finished. 


46. We now discuss the convergence situation in an example. Consider the 
matrix 


(177) A=(_3 oye 


which has the eigenvalues 1 and 2. The eigenvectors corresponding to the eigen- 
value 2 are £=(1, —1) and 7 = (1, 2). 
We shall discuss the convergence of A, to 2. Put 


(178) eats oes & = 6 (45%) 5 A= 11, 9), 


(179) %,= —1+49,, Vn = 3 — On: 


Then we have 


4— % Vx t2%y,—3V x 
BE ue) = 144, Vy 7 


4+ (1—o,) (§— ex) —2(1—9,) —3 (3—@x) 


— 24+46,460x+3 Fx On 
1+26,4+30.—36, 0x” 


(180) [pk 9x Rx 


Further, from (78), using (178) we have 


es 3+(3+Mx) Fet1 
es 2 by + 2% y44 
or 
; a slap Emel) Xx 
Xy 41 


—(3+4)—24, 


iS) 
Go 
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and, using (179), 
3= (2H) (1 = 4) 4 


Oy = 3 tie 2(1 G,) : 
(181) ~ Oe 
Pit yt 20, 
In the same way, 
y an 2—(3+u,) Vut+1 
@ 2 3 Veta ‘ 
*TE (3+ by) +3 Vx’ 
and again using (179), 
. als Ox Ux ’ 
(182) Ox+1 ian Ne ares 3 Ow ‘ 


47. If we start-with £,=7,;=(1,.1)¢,=2, @o=— 4> from (180) weshave 
=6 
= — =1; here at the first step we get an eigenvalue, 
Morea eran eS 1, Ao r p we g g 
although &) and mp are not eigenvectors at all. 


If we put 
Ej == (1,0), %=V=0, =1, Oo= 35> 
we obtain from (480), (181) and (182) the values of o,, 0,, M,: 


Oy = .07582938; 05 = .050552923; lg = .026707002, 
3 = .0017186302; 03 = .0011457538; fg = .0417601246, 
6, = .0°30145883;  04= .072009726; ug = -0"45.452646. 


If we compute the quotients 


Pichit 8 129 Page 


(A,—0)3 be 
for x =1, 2,3 we obtain the following values taken to 4 significant figures: 
35793 92373" 41: 


If we now apply to the sequence A, the argument in Section 19 Part II of this 
paper, we can expect an improvement upon A, and wv, if the formula (61) is used 
with x =3. Thus we find that 


Ay Jp = = 01476, et See Potton a73 
A, al An ae .026707 eeey (Ay =; Ash —— ha .0419049 eee 
ji, fy Ses As Clacce eM eee 


.0419049 
which is better by one decimal than A4,=2.0'4545 .... 
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On the Convergence 
of the Rayleigh Quotient Iteration for the Computation 


of the Characteristic Roots and Vectors. IV 
(Generalized Rayleigh Quotient for Nonlinear Elementary Divisors ) 


A. M. OSTROWSKI 


48. In this part * we shall give a quadratically convergent iteration rule for 
computing an eigenvalue of a matrix to which there corresponds a non-linear 
elementary divisor. Our method generalizes the rule discussed in Sections 1—8, 
Part I, to non-symmetric matrices and combines it with rules for hastening the 
convergence of an iteration of the first order. 

Let A be an Xm matrix with an eigenvalue o, and assume that the maximal 
exponent L of an elementary divisor of A corresponding to o, is >1. Choose 
two (row) vectors, « and f, and consider for a A) which is different from all eigen- 
values of A the linear equations 


(183) (AAI) =a, n(A— Al) =B, 
defining a couple of vectors &, 7. 
Now form the generalized Rayleigh quotient of € and 7 as defined by (72), 
and put 
(184) Ay = R (Em) =, 
assuming that 
(184°) n&' == 0. 


Then from (184) and (183) we obtain A, as a rational function of Ay, 


(185) A, = 9 (Ao): 


We shall have to prove first that, under suitable hypotheses, the iteration 
by the function p(A) defined by (185) has o as a point of attraction. This is a 
consequence of the following 


Theorem. In the notation and under the assumptions of Section 48, there exists 
a matrix H depending only on A and on the choice of o, such that if « and B are 
 * This paper was prepared in part at Numerical Analysis Research, University of 
California at Los Angeles, and in part under a National Bureau of Standards contract 


with The American University, Washington, D.C., with the sponsorship of the Office 
of Naval Research. I am grateful to Mr. G. CULLER and Dr. E. V. HaynswortH 


f or discussions. 
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a 


chosen to satisfy 


(186) pHa’ +0, 
then we have 

, 1 
(187) y'(o)=1—-,. 


In particular, if (186) is satisfied, then the condition (184°) 1s also satisfied when 
A, is sufficiently near to o. 

49. (187) shows that the iteration by the function g(A) converges rather 
slowly, though linearly; even for L =2 the derivative is $. However, since p (A) 
is a rational function, it follows from (187) ae as soon as (186) is satisfied, we 
have a development 


(188) ei) =o+(1—7 J (A ~a) +> (= of, 


and therefore different methods of acceleration are applicable. 
Consider for an «+1 the function 


(189) p*(A) == (pa) — 4). 


We verify at once that y*(a) =o and obtain for the value of the derivative 
of g* ing: 


(190) Fl ens 


If now the value of L is known, we can choose « =1 —+ and obtain for the 
corresponding function @*: 


(194) pi(4) =L (a) —(L—1) 4 
with vanishing derivative at the point o. The iteration by means of this function 
then converges quadratically to the value o. 

If the value of L is not known, then it is still best to use the iterating function 


(192) (A) = 2(A) — A. 


The iteration by this function converges quadratically for L =2, while in the 
case of a general L we obtain a derivative 


(193) 73 (0) et 


which is less than 1-+=9'(6 ). On the other hand, if L>2, the value of L 


will usually be recognized easily after a certain number of steps, and then the 
corresponding function , (A) can be formed. 

The use of g.(A) is best under the assumption that the value L =2 is the 
most probable of all L>1. On the other hand, if all values of L between 2 and 
m are more or less equally probable, we shall take 


1-45 
ie : 
im 


(194) 
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and obtain readily 


(195) ky ieee 
i 


50. An iteration with quadratic convergence can be obtained in our case, if 
L is not known, by the Steffensen method, which, however, implies a double 
amount of computational work, for it uses a combination of 2, p(A) and p(p(A)). 
In this method we form the function 


_— — o (p(A))— 9 (A)? 
(196) P(A) A—29(A) + (9(A) ’ 


which for A= A, usually gives a better approximation than 
Ay = 9 (Ay) = P ( (Ao) - 


In practice the value of y(d,) is of course computed by obtaining successively 
the values of A, and A, and by taking 


s yi Ree 
Cy P(A0) Py — 2A thy 
If then we put 
(198) Ay = P(A), 40 = Ap, 


we can consider the passage from jy to Ay as a direct iteration by the iterating 
function ®(A). 
However, this iteration converges quadratically. We prove this by using the result 


of our paper*. If we replace there A, and A, by 2, yby A, Cbyoanda,=a, by 1 —; : 
we obtain from formula (21) /.c. 


PNAC Mag ZA 
aoe = L?T,+0(A— 0) 
where 7, is obtained from the formulae (15) and (9) J.c.: 
4 N2 il ' 
a (1 = z) E (p(a)) — (1 i a E(A), 


E(A) = a, + 43(4 —0) +, 


(199) 8) 5 ARG ig O (0): 


The one theoretically unsurmountable difficulty in the practical application 
of this method appears to be the fact that the bilinear form in (186) is unknown 


* Ostrowsk!, A.: Uber Verfahren von STEFFENSEN und HovusEHOLDER zur 
Konvergenzverbesserung von Iterationen [Mauro PiconeE zum 70. Geburtstag, 
ZAMP Vol. VII, 218—229 (1956)]. 
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so long as the transformation of A to the Jordan canonical form has not been 
carried out, and this presupposes the knowledge of the eigenvalues. However, 
in practice this is hardly a difficulty at all, since it is infinitely improbable that 
for a and B taken at random, (186) would not be satisfied. 


51. Lemma 1. Consider the matrix 
(200) A,=ol+U, 


I being the unity matrix of order | and U the corresponding auxiliary unity matrix 
of order 1, which has 1’s in the first swperdiagonal and zeros elsewhere. Then we 
have, for Ao, A>c: 


e he ae 1 
(201) (A,—Al) etter Oa) 

tee Ee, 1 
(202) Udy At) t= td Opa ee 


Proof. Putting x = cs we have, since U'=0, 


(6 =A) I-40)? =I — x0)? 


f—=1 
= aoe ye 7 Us 
y=0 
and, multiplying this by U, 
1—2 


U((o — A) I+U)2# => (v +1) 2? UH, 


v=0 
as U'=0. Taking out the highest terms on the right, we obtain (201) and (202). 
52. Lemma 2. Let o be an eigenvalue of the matrix A to which correspond 


elementary divisors with the maximal exponent L>1. Then there exists a matrix 
H =0 such that we have for A«: 


‘y H 1 
(203) (A A1) z= L oee eae 
29 Tye Seed Ate Ae ! 
(204) A(A—AD?-—o(A—AN?=(L—1) Teer aah ar esse 


53. Proof. Since the assertions (203) and (204) are invariant with respect 
to a similarity transformation, though H does change, we can assume from the 


beginning that A is given in the Jordan canonical form. We can therefore write 
A as a direct sum, 


(205) A=)°A,+B, 


where B is a matrix with eigenvalues =o, while each A; is a matrix of order m, 
given by 


(206) A, = Lat U,,, (1 =1,...,k), L=maxm,, 


7 


I,,, and U,,, having meanings analogous to those of J and U in Lemma 1. 
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54. Then for A->o we obviously have 


(Ben ot O'(1) 
and therefore 
k 
(eA) Ae). ee OA). 
i=l 
Applying the result of Lemma 1 to each of the first k terms in the direct sum 
on the right, we obtain 
ae Te 


5 eae i : ersly PS 
i i bs ae ees: dhgenak 


Here in the direct sum 
(207) =a 


‘ : L= 
we must sum over all 7 with m,;=L. However, for m,<L we have U,, 1—0, 
and therefore we can write 


where H +0, since there exist m,=L. 
This proves (203). 


55. Further, from (205) we have 
h 
A(A —AD)*—o(A —AN*= > (A; —o Ip) (A — AT y+ OU). 
=A 


Applying (202) to the first k terms in the sum on the right, we have now 


Um -1 4 
= ——— 3 4 Ae 72 —— yy. ihe | O = ; 
(A; Gils (A; A pss) Ura (A; AI;) (m; 1) (A—o)™ ates 
Therefore 
DU ee 
i 1 
a -2 =2 | : A 
Ni 2 6A A 1) aE ol reece 


where the numerator of the first fraction on the right is given by (207). This 
proves (204). 

56. Proof of the theorem of Section 48. Take the matrix H occurring in 
the formulae (203), (204) of Lemma 2 and assume that the condition (186) is 
satisfied. Then from (184) and (185) for Ayo we have 

ete gee 


From (203) and (204) 


ie = L—1+0(A,—<o) 
(oe *t © treed) 


346 A. M. OsTROWSKI: 


This proves assertion (187) of the theorem. The last part of the theorem follows 
immediately by (203) from 
Pap te Se 
nk’ =B(A — Ag) *a eater tO geeer clk 


57. We consider as an example the matrix 


Zibwat 
1-( 1 s) 
{4 


Here (A —AI) has the elementary divisors (A—1)? and (A—2), and we have 


D=|AS All = AS1Y CAA aA SA 2) 
We easily obtain 
A2— 21+ 4, =a 324-3 
D(A~AIA= Vee Via Ane eee 
A— 2, 2—A, 22—3A442 
.. HOE es = 
cay marae Ane? 
AB — 342 +124 — 16, =O) 245" 6727 —15 2-9 
-( 2/7 + 44 —12, AB—5A?—A+11, 62? —15449 
222 — 84+ 8, —2A7+ 8A—8, A—527274+ 81-4, 
and 
i aa A ae 
2 (A? — 4), — 6(A2 — A — 1), 3 (As — A? — 20 + 2) ' 
=| 424+ 5/7—8A—4, A112? +14142, Grit | 
AB — 322 + 4, —(A— 322+ 4), A®—5 22481 —4) 
From the definition in Section 48, 
2 — BYa’ 
(208) p(t) =F. 


Since we have 


A® shtidcni 8 eS (Qik BN Sidi ele 
yoxaaon{ wane — 3 (24 — 3) sa" |, 
(4 — 2)? ed ey 0) 
it follows that 
(209) O(A)—1 _ BZa’ 


At BXa’’ 


For A=1, Z and X become respectively 


gar oman? | 10) == 6 \n-6_ a0 
Z = aS 3 O ) Xy= (5 6 0) = 22, 
i = 7h 30) =7 0 
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and we see that as A->1 the ratio (209) tends to 4 if BZ,)«’+-0. Z, is therefore 
the matrix H in the condition (186), and if we write 


& = (A,,4g,43), B= (dy, dg, bs), 
this condition becomes 


(4, — 2) (— 3b, — 3.2 + bs) += 0. 
We see that « and # have to be chosen in such a way that a,=-a,, b3=+=3(b, +05). 


If, for instance, we take 


e050), p= (1,0, 0): 
we have from (208) 


VA 
soe Cant, 
(210) vA) Aye 
and from (192) 
waite Herat 
(241) pr(a) = AE 
If we start with 2,=0 we obtain from (210) the following values of 2, and ; s foe 
besa! 
glenda 
4 ‘ 1—A,-4 
6) 0. 
4 4 6 
2 6526 579 
3 .80648 557 
4 89599 5375 
5 94565 5225 
6 97214 5126 
Starting again with 2)=0 we obtain from (211) 
(1 i Ay—1)” 
A, 1—A, 
O 
8 5 


.984616 2.046 
.999884338 2.03347 
.98331229 2.072006 


AOwNwW |} CO ®@ 


Certenago-Montagnola 
Lugano, Svizzera 


(Received February 2, 1959) 


Characterizations of Laplace Transforms 
and Perfect Operators 


JaDAWESION 


Communicated by A. ERDELYI 


1. The notion of a perfect operator was introduced in [/], mainly for the 
purpose of widening the Laplace-transform interpretation of HEAVISIDE’S symbolic 
calculus. A problem mentioned in [1], and related to the classical problem of 
characterizing Laplace transforms, is that of finding a necessary and sufficient 
condition for an analytic function to be the Laplace transform of a perfect operator. 
We give here a solution of this problem, in a way that enables us to show (in 
answer to a question raised in [2]) that all perfect operators are obtainable by 
repeated “‘differentiation”’ of continuous functions. We also discuss some impli- 
cations of these results. The notation and terminology are those of [2], supple- 
mented where necessary. 


2. We shall require a pair of lemmata concerning Laplace transforms of 
functions: we give proofs of these for the sake of completeness, though in fact 
much stronger results are known. 


Lemma 1. Let x be a continuous function of a positive real variable, and suppose 
that, for some real number a, e~“'x(t) is bounded when t>0. Let x be the Laplace 
transform of x, and let c be a real number greater than a. Then x(z) is bounded 
when Rz=c. 


Proof. When & zc, 


so that, if |e~*’« (t)| Sk when ¢>0, 


FS fle“ xl dts f kee-o'gy—_* 
0 0 


Cc—a’ 


Lemma 2. Let g be a function which is holomorphic in a domain containing 
a closed half-plane Rz=c, and suppose that 22g(z) is bounded when Rz=c. Then 
the formula Ns m 

ect f 4 
«Q)= 5 [ gle+ia) MSE (¢ > 0) 
—0o 

defines a continuous function x such that e~°'x(t) is bounded when t>0, and 
% (2) =¢(z) when Bz>c. ; 
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Proof. Since (c-++-1w)* g(c-+im) is bounded when @ is real, 


x(t)|S eiifle (c+iw)|dw< o; 


and there is no difficulty in showing that x is a continuous function. If #z>c, 
then 


(2) = > de (om "fe (ct+io)é°'dwdt 


[o.e) 


=f tio) [« eletio—2)t digg — 1 Ae) dw 
2% J Chiw—z . 


—co 


Now, for each positive Sent e, let S, and C, be, respectively, the straight 
path and the right- ee, semi-circular path from c—7zo toc+7o. By the residue 


amigt) =f £8) as— f #8 a 


Se 


But, when s is on C,, sg(s) is uniformly bounded with respect to 9; so that, 


by a standard argument, 


Thus, when #z>c, 


g(z) = lim ! AOS Geen! i dw =X (z) 
oc 2ni J s—z 2x J ction—z aan 
So —oo 


3. We recall that a perfect function, of a positive real variable, is a function p 
which has derivatives of all orders and is such that, for »=0, 1, 2,..., p™ (2) 
is of exponential order (as too) and p”(0+)=0. The following result gives 
a characterization of the Laplace transforms of perfect functions. 


Theorem 1. A function p, holomorphic in a half-plane Az>a, ts the Laplace 
transform of a perfect function if and only if, for every integer n, there 1s a real 
number c, such that 2" p(z) is bounded when Rz=c,,. * 

Proof. Suppose first that qm is the Laplace transform of a perfect function f, 
and let ~ be any non-negative integer. If x=p” then ¥(z)=2z”"q@(z), so that, 
by Lemma 1, there is a real number c,, such that 2” (z) is bounded when #z=c,,. 
We may take c,, to be positive, so that z~”q(z) is also bounded when #z2c,. 

Conversely, suppose that, for every integer 7, there is a c, such that 2” (z) 
is bounded when #z=c,. Then 2”*2@(z) is bounded when #z=c,,,,, so that, 
by Lemma 2 (with g(z)=2"(z)), there is a continuous function x,,, of exponential 
order, such that pad 
Xn (2) = 2" p(2). 

* It is not enough that 2”f(z) should be bounded when #z=c,, for every m: this 
condition is satisfied if f(z) =I'(z), for example, but then 2”/(z) is unbounded when 
Az=c, for every real c. 
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Now, since %, (2) =272 X41 (4), 
t 
XQ) =JS %nii(u)du (> 0), 
0 


so that x,(0+-)=0 and x,,,=%,. Thus % is a perfect function, with x) — x, 
when 0. Since y=Xo, the theorem follows. 


4. Our main result is 


Theorem 2. A function f, holomorphic in a half-plane Az>a, 1s the Laplace 
transform of a perfect operator if and only if, for some integer n, there is a real 
number c such that 2"f(z) is bounded when Az=c. Every perfect operator ts of the 
form D”™X, where D is the differential operator (on perfect functions), m ts an 
integer, and X is the convolution operator determined by a continuous function. 

Proof. Suppose first that / is the Laplace transform of a perfect operator f(D), 
and let us assume that, for every integer and every real number c, 2”f(z) is 


unbounded when #z=c. Then a sequence 2, 22, 23, ... can be defined inductively 
soithat:Z2,—-¢ and, for w= 14245) ue, 

(1) len” #(z,)| 2m, 

(11) Rigi = Callan 


Let w,= log z,, for each m (the logarithms having their principal values). Then 
it follows from (ii) that |w,,,/w,| 22, and hence that the infinite products 


[e.°) 


oo 

Pape —2w/|wp| 

Tes og eae 
n=1 


n=1 


converge for all values of the complex variable w, and represent entire functions. 
For each positive integer m, let y,, be the entire function defined by the formula 


zs pene Su | | w —2v/|w,| 
ie (w) (1 Wy (1 = i: q 
Then 


(iii) Ym (w,) =O (m = n) 


and, for each m, 
m—1 0O Paice 
vse (FF S|} 0-H 


CO 
where ee Nall Now |w,,/w,|22 when nXm—1, and | w,,/@,|2"-" 
n= 


when »=m-+-1; hence 


lynn 2{ TP (1—20—™)| Jemeem], 
n=m+1 


Thus 


(iv) | Wm (7e,,,)| = ) | 2m ae 


lo) 
where ae (1—2°"), a positive number which, like a, is independent of m. 
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With 2z2e, let w= log z. Then Zw=1 and | Jw|<2/2; so that 2Zw=|w|. 
Hence, for each 2, 


(2 — 2) geet 


Wy 


s(i+|z 


Miles s (1 as Pall e7lelen| = 4. 
Wy 


Since, for each m, 


Wp (0) = eW2#lleml TT (1 cat =] e~2lltnl 


n+=mM On 


it follows that, when #z—e, 
(v) [vn (w)| S14. 


INOW, 10% 77—4, 25 3, . 2. and Be 2c) let 


Pm (2) = eV Oz) y,, (log z) . 


Then each function @,, is holomorphic in the half-plane #z>e, and, for each 
non-negative integer 7, 


|" gm (2)| S| zl" | eV] | py (0)| fa] em H22P90, 


by. (v)., But, if r= |/|'2/z,, 


, 


| z|” | ea) allel) <& | Z|” qn eo? 


’ 


and 1?”e~", as a function of the positive variable t, has the maximum value 
(2n/e)?”, attained when t=2n. Thus 


(vi) | 2" Gm (2)| S10?" | Z|": 


It is a consequence of (ii) that, for any integer n, 


lee) 


(vii) a Lemitcna < oO. 


m=1 


Let 


p (2) = 2 Fg hil OC? ot 


By (vi) and (vii), with »=0, this series converges uniformly when #zZe; so 
that the function y is holomorphic in the half-plane #z>e. Moreover, when 
n=O, 
lo, 2) 
|2e"p(z)| Sn" Bile ele 

Thus, for every integer 1, 2” (z) is bounded when #z2e. By Theorem 1, there- 
fore, p is the Laplace transform of a perfect function, say #. 

Let q be the perfect function /(D)p. Then, when #z is sufficiently large, 
9 (z)=f(z)@(z). In particular, for each positive integer m, 


q (Zn) a AC) | a AP aie Pm (Zn) , 


352 J. D. WEsTON: 


since Pm»(z,)=0 when mtn, by (iii). Hence, by (i) and (iv), 
1F (Zu) | S 2 | Zu |? | enV ml) | | Z| 2S 10 OV? 6. 


Thus 7 (z,,) > co as moo. By (ii), therefore, 7 (z) is unbounded when #z=c, 
for any real number c. This contradicts Theorem 1, forcing us to conclude that, 
contrary to our initial assumption, z”/(z) is bounded when ARz=c, for some 
integer 2 and some real number c. 

Conversely, suppose / to be such that 2”f(z) is bounded when #z2c. By 
Lemma 2 (with g(z)=2" *f(z)), there is a continuous function x, of exponential 
order, such that 


(2) =e 2) 


Thus / is the Laplace transform of the perfect operator D” X, where m=2—n 
and X is the convolution operator determined by x. This completes the proof 
of Theorem 2. 


5. In [2] it is shown that the perfect operators can be characterized as those 
linear operators, on perfect functions, which commute with the differential 
operator D and are continuous in a certain sense. We take this occasion to point 
out that not all linear operators that commute with D are perfect. To this end, 
we note that the perfect operators whose Laplace transforms are rational functions 
constitute a field, %, and that the perfect functions constitute a vector space 
over 7, in an obvious way. This vector space is not one-dimensional, and therefore 
has a basis containing at least two perfect functions, p,q; and we can define 
a linear operator A on the space in such a way that df=0 and A4qg=+0. Such 
an operator cannot be perfect, since the algebra of Laplace transforms has no 
divisors of zero. However, since A is linear with respect to %, it commutes with 
D (which belongs to %) and is linear with respect to the complex numbers (which 
are embedded in %, in the sense that %} contains the scalar multiples of the unit 
operator, J). 


6. Systems of linear differential equations that occur in certain physical 
problems can be formulated in terms of perfect operators: all functions of the 
time variable are replaced by the corresponding perfect operators (which may 
then be taken to represent “generalized functions’’, including impulses if neces- 
sary), the symbol “dx/dt’’ being replaced by “DX —x(0+)I’, and so forth. 
If a solution can then be found within the system of perfect operators, this will 
have a direct physical interpretation, since, by Theorem 2, every perfect operator 
represents an ‘‘m' derivative’ of a continuous function. Thus “spurious” 
solutions cannot occur. In the system of operational calculus which has been 
developed by MixusrXsxkr [3], solutions sometimes occur which have to be 
rejected on grounds that seem extraneous to the original mathematical frame- 
work. For example, in the discussion of a certain type of boundary-value problem, 
the expressions “el?” and “‘e~!”” appear to have equal validity, although only 
the second has a physical meaning. In the system of perfect operators, on the 
other hand, it is only the second expression that has a mathematical meaning, 
since the function el? does not satisfy the conditions of Theorem 2 (in fact it 
was shown in [/] that this function is not a Laplace transform) 
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In MIKUSINSKI’s system, the “operators” constitute a field; but it is easy 
to see that for some purposes a field is too large. For instance, the operator e~“”, 
where «> 0, represents a “‘pure delay’’, so that, since there are functions which 
are not obtainable by delaying others, it is inappropriate that e~“? should have 
a reciprocal in the system. It is of course possible to restrict attention to those 
elements of the field which represent generalized derivatives of continuous 

functions. In the system of perfect operators, this restriction is inherent; and 

there is also the advantage that each continuous function considered has a 
Laplace transform, an “operational form’’ of the function being obtainable 
by simply replacing the symbol for the complex variable, in an expression of 
the Laplace transform, by the symbol “D’’. Moreover, Theorem 2 provides a 
simple criterion for deciding whether a given analytic function is a Laplace 
transform or not; and if it is, we can “invert”’ it by ‘differentiating’, a suitable 
number of times, a continuous function which is given by a Fourier integral 
according to Lemma 2. The price of these advantages, namely the condition 
“x (t) = O(e°’) as t—> co’, is not exorbitant. 


The field %, consisting of all perfect operators that can be expressed rationally 
in terms of D, contains the solutions of all finite systems of ordinary linear diffe- 
rential equations in which the coefficients on the left-hand side are constant 
and the “functions” on the right-hand side have rational Laplace transforms. 
After formulating the equations in terms of perfect operators, the solutions can 
be found by the usual processes of rational algebra (which are valid for any field), 
and interpreted in terms of partial fractions. 


7. The use of perfect operators in linear electric-circuit theory is closely 
analogous to the use of complex numbers to represent currents, voltages, impe- 
dances, and so on, in the steady-state sinusoidal theory with which engineers 
are familiar. The various “operational impedances” which determine the response 
of a network to any given excitation can be found as follows: first calculate 
the steady-state impedances in terms of a frequency parameter y, and then sub- 
stitute “D/2a7” for “vy” at each occurrence. The operational impedances of a 
finite network are always elements of the field %, and all the usual network 
theorems remain valid. 


Differences of opinion have been sharply expressed, even in quite recent 
times, over the relative merits of the “operational’’ and the ‘“‘transform’’ points 
of view in interpreting HEAVISIDE’s symbolism (see, for example, [4] and [4]). 
It may be hoped that these differences can be reconciled by appealing to the 
natural isomorphism (denoted in [2] by “’’) between the algebra of perfect 
operators and the algebra of Laplace transforms. However, in the notation used 
by HEAvISIDE, and still sometimes advocated, the “operational form” of the 
unit function is “1’’: this is unnatural from the operational point of view that 
we have been considering, since the perfect operator determined by the unit 
function is the integrator, D+, not the unit operator J (which represents a unit 
impulse at the instant ¢=0). Our terminology could be altered to accommodate 
HEAVISIDE’s notation, with perhaps some practical advantage; but then the 
operators “determined” by functions would not form a ring, and their Laplace 
transforms would differ (by the factor z) from those of the functions themselves. 
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The Cauchy Problem and the Mixed Boundary 
Value Problem for a Non-Linear Hyperbolic Partial 
Differential Equation in Two Independent Variables 


JAMES CONLAN 


Communicated by L. CESARI 


Preface 


The following paper deals with the Cauchy problem and the mixed value 
problem for a canonical form of a second order partial differential equation of 
hyperbolic type in two independent variables. The exact statements of the 
Cauchy problem and of the mixed value problem dealt with are given in Theo- 
rems 1 and 8 respectively. 

Theorem 1, which concerns the Cauchy problem, has been proved by 
P. LEEHEY [1], see also P. HARTMAN & A. WINTNER [2]. These two theorems are 
proved here by a variation of a method employed by J. B. Diaz (to whom the 
author is indebted for his assistance in the preparation of this paper) in the 
treatment of the characteristic boundary value problem [3]. The main interest 
of the present approach is that the existence theorems are proved in a way 
which leads directly to a method for the actual numerical construction of solutions. 

Appendix 1 is concerned with concrete numerical applications of the present 
method to the mixed boundary value problem and to the Cauchy problem. 
In order to make the paper self-contained, Appendix 2 presents a slight generaliza- 
tion of a relatively little known theorem of ARZELA, which is fundamental for 
the present method. 
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Section I. The Cauchy Problem 
1. Statement of the Problem 
Section I of this paper deals with the Cauchy problem for a hyperbolic partial 
differential equation in two independent variables, the precise statement of the 
problem being given in the following theorem. 


Theorem 1. //: 

a) fis areal valued function of the five real variables x, y, 2, p,q (this includes 
the case where f is independent of some or all of the variables x, y, 2, ~, gq); 

b) f ts continuous for OSx, ySA; —w<2,f,9<00; 

c) There is a constant M=O such that |f(x, y, 2, ~,q)|SM for 0Sx, ySA; 
— 0 <2,f,9< 00; . 

d) f satisfies a Lipschitz condition in p and q, 1.e., there 1s a constant L>O 
such that for OS %, ySA and —~w<z, f, q, by, I< we have |f(x, v, z, , g) — 
f(%, ¥, 2 Po» W)| SL{|P — do| | = Gale. 


e) o and t are continuous real valued functions of the real variable 2 for OS ASA; 


then there 1s at least one u such that for OS x, ySA: 

f) wis a function of x and y; 

g) u and its partial derivatives u,,U, and uU,,=Uy, are continuous ; 

h) tyy(%*,¥) =f1%, ¥, U(%, Y), Uy(%, VY), My(%, 9); 

1) u(0, 0) = constant (which can be taken to be zero); 

j) U,(%, ¥)|yax = (4); 

h) u,(%, ¥)|yax= T(%). 

The hypotheses of Theorem 1 are not sufficient to guarantee the uniqueness 
of the solution, u(x, y), as the following example shows. Let f(x, y, u, p, 9) 
=—|u|*, u,,=—|u|#, 40, 0)=0, o(x) =0, t(x)=0. Then u(x, y)=0 and 
u(x, ¥) = yaa (vy — x)fare both solutions. However, it is known that if 1 (%,9,2,P59) 
satisfies a Lipschitz condition in z, in addition to satisfying the other hypotheses 
of Theorem 1, then the solution, u(x, y), is unique. This is proved in [5] for the 
linear case, but the proof is the same in the non-linear case. 


2. The Approximate Solution 
If f is a function of x, y alone, and is independent of z, f, g, then a solution 
of CAuCHY’s problem is given by 


(2.1) Hoy) = fay fHEn) ae +] ol8) a8 +f 218) a8. 


That (2.1) is a solution is verified as follows: Clearly “u(0,0) =0. Also 
y 


u,(x,¥) =J f(x, n) dy +o(x), 


and so “,(x, x) =o(x). Similarly 


Uy (x,y) = SHE: y) d&+1(y); 


hence #,(%, x) =t(x). Finally 


(2.2) Uyy (XY) = Uy (%, ¥) =F (x,y). 
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For a given positive integer 1, let the square R (OX x, yS A) be subdivided 
into n* rectangles M7, (7 =0,...,%—1; k=0,...,~—1) as follows: Choose n+1 
real numbers. £; (7 =0,,..., »)s0ch that 0O= <<. <f ,.<i/=A. Then the 
vertices of 37, will be (¢7, th), (41, te), (GF, Hai), and (41, 41). 

We shall use the following notation: 


We ={(% VG << Gas <y < hiya}, 


Meee) ob See SV hs}, 


VF =} (4; 9) 1G Sesh; BSyStai}, 


An 


etc. We shall also let (x,;, y,) stand for a fixed but arbitrary point in Yi ,. 

It can be readily verified that formula (2.1) would still give “‘solutions’”’ to 
the Cauchy problem for (x, y) in the interior of Kj, if f(x, y) is allowed to be 
piece-wise constant on }t (7.e., to assume a suitable constant value in each square 
Jt?,, and thus to allow the possibility of discontinuities along the lines x =#,; 
and y=+#,); while o and 7 are replaced by piece-wise constants, 7.e., constant 
for FSFS 4,44, 1=0,-.9, t= t--This fact will now be used to derive a scheme 
for the approximation to the solutions of the Cauchy problem which will be of 
basic importance in proving Theorem 1. 


Our “approximate solution”’ is then 


y, U] x Ed 
(2.3) u,,(x,y) = —J dn J F,(E,n) 46 val 2, (§) dé ral T,,(§) dé, 


where the F,, 2',, and 7, are defined by the following expressions (where we 
have dropped the subscript and superscript ”, as we shall continue to do whenever 


no confusion will result): | _ o(t,); t=t(t,) 
7 ql? ) q/* 


T(x) es Opty if be St, 415 
ie 09 if Xx =e 
T= 1G) if be % <heays 
4 Ty if x=t,. 
[tall t,) if 7+ k; 
UW; p= \ 4 4 aes 
| fEdE+s Tae if F=h. 
0 0 
[tiene if k>7 and (%, y) ER, a; 
v(x, ¥) =) %, j iA YEN 7: eS 
exe ifk<j and (x,y) €|Rj,- 
es ifk>j and (x, y/ER, a]; 
2(x,y) = fn if (x, y)ER; 53 sae 
t, ifk<j and (x,y)€|R;,- 
» if k>j7 and (x,y)ER, x]; 
W(X, 9) = ; (+441) if (x, 9) EN, 5; mike 
ist ifk<j and («,y)€[Rj,- 
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veh if k>7+1 and (x,y) EM, al 
teubis— Mase if h<j—4 and (x,y) €[B a) 
bry 
0 if (x,y) EN, 53 
0; if (xe Ry_1 | UR, ja 
tithe 1,k if k>7+1 and (x, ¥) ER, als 
R+1 a 
i | nett ifk<j—1 and (x,y)€[Bi,; 
tr; if t( y) ER 1, | UPR, ja 
0) if (4; VM Fs 
f[A(x, ¥), w(x 9), 9 (x, 9), (x, 9), 0(% ¥)] 


At first glance our definitions appear to be circular, since w,,(x, y) is defined 
in terms of F which is itself a function of the w, ,’s, which in turn were defined 


in terms of w,,. 


However, it is seen that starting from the line x =y the “‘ap- 


proximate solution”’ (2.3) can be continued into successive }; ,’s in a well defined 
manner, 7.e., such that only previously computed expressions occur as arguments 
of F, and hence is well defined in all of #t. To show this explicitly consider the 


following : 


If (x,y) CR, ;, then by (2.3) 


=f 2 Natt fT a 


If7<k and (x, y) CR | then 


u(x, ¥) = 


= fi, Tele tua) (y 


— fan J FEn) ae +f BE Vac ef Menge 
0 


y veal 


= — fd J FE.) dé — 


— Jay [FG m) a8 +f 28) a8 +J T(E) dé — 


ba 
~fan IPE, ) )dé + J X16) ae +f TU) é) dé + 
Ui 1) vee 


+ fan J Fn) dé —J 6 def TI@) 4 


Gir tty 


— t) + 4(%, t) + U(ti41, y) = Uj 4 1,k> 


and hence u(x, y) is determined in j;,, in terms of functions already computed 
in adjoining rectangles. 
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3. Some Properties of u,, 


From (2.3) it is seen that if (x,y) ER, 


OU . 
an (x, 9) = NO Ga ees 


(3.1) 


OUn 


ay ee ~ fPG,9) 48 + Ty). 


By the definition of F, XY, T, it is seen that if (*, V) ER,» 


Oru Ou, Cu 
ae (t;, y) = ta (x, ¥) Se (t;41,), 
(3.2) 
On, ( t,) Sr OUn (x y) oO Uy, (x, ¢ 
ay a CMe a OY »"e14) > 
where, as usual, 
Ot Uy Un(*%+h, y) —Uy(%, ¥) 
XG sli g NAB SA) 
Ox ( »¥) pay x h ? 
and 
oun = 14 Un (~—h, y) Uy (%, y) 
Ieee a0 |, al =F 3 


Another useful relation is 


Muh Toh sty [ fan fae—fen fal 76 


titi 


41 ta 


va fie faae 


aa 


ty tha 
5 J an [PEM dE to 


es hj 


tk 
= f F(x.) dé +oj41 


baa 


(3.3) wig oe: 


Oo Uy OUn OU, 
= ; (41) ¢a) ao Ox (x;, t) = (¢;, t,) . 


Ox 


We let «, = max |¢? —t¢*,,|. We are going to consider a sequence of sub- 


jsn-1 
divisions of it oxen that «,>0 as > o, and we wish to show that the sequence 


{u,,+ contains a subsequence which converges to a solution of CAucHy’s problem 
As a first step we have the following two theorems. 


Theorem 2. The elements of the sequence {u,} are uniformly bounded. 


Proof. Since |/(%, , 2, ~,9¢)|SM, the same inequality holds for F(x, y). 
Similarly since o(x), t(%) are continuous on 0<%SX4, they are also bounded 
in absolute value there, and hence so are (x) and T(x); we can without loss of 
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generality assume these bounds to be M. Thus for 0S%, ySA we have 
F(x, <M, [Z)/SM, [Td] SM. 


Hence by (2.3) |w,(%, y)| S(4?+2A4) M, q.e.d. 
Theorem 3. The sequence {u,} is an equicontinuous family of functions in the 
region OS %, YSA. 


Proof. We shall consider the case where x< x, y<¥y. The treatment of the 
other three cases (namely *+<%, y>y; x>x, y<y; and x>%, y>y) will be 
clear once this case is disposed of. 


[ns (F5) — males 91S { fdr, fads— fan fash [Fem] +f BOA +S (TE) a8 


cf Pdejons fafa nents fineie sini 


max(y, #) 


<M{A|x—z|+4|y—9| +12 


q7e.d; 


Now Theorems 2 and 3 together with AscoLi’s theorem (see Appendix 2) 
imply that there is a subsequence {u, } of {w,} which converges uniformly to a 
continuous function u in the region X. 

In the special case that f is a function of x, y, z only, but not of # or gq, it 
would be a straightforward task to show that w satisfies an integral equation 
of the form (2.1) but with /(é, 7) replaced by /[é, 7, u(&, 7) |. However, instead 
of carrying out the details of this proof, we shall proceed to the general case 
where / is a function of all five variables x, y, z, p, q. 

The added complication in this more general case is that while the two first 
partial derivatives of w, exist in the interior of each region Rt; , (OSj, RSn—1), 
the partial with respect to x of w,, need not exist on any of the lines x =¢,, and 
similarly the partial with respect to y of w, need not exist on the lines y =1;. 
Since these lines become dense in the region {t as 7—> co, ASCOLI’s theorem gives 
us no information about the existence of convergent subsequences of {@u,,/0x} 
and {éu,/dy}. 


Before proceeding further we note that by (3.1) 


ie (x,y )]sM {ly —2| +1}<SM{A +1}, 


and similarly for |du,/@y|. Hence a(x, y) and e(x, y) are uniformly bounded 
for (x, y) ER. 
4. A Useful Inequality 


The following theorem will be of fundamental importance in circumventing 
the above mentioned difficulty. 


Theorem 4. /}: 
a) tis a non-negative integer ; 
D) 80s Bis -++) & IS a sequence of non-negative real numbers: 


C) Rae e. span tsa non-decreasing sequence of veal numbers; 
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k—-1 


ad) there are numbers L>0O and e>0 such that g,Se+L > g,(2,.4—2,) for 
k—1=r, r-+1,...,t; then re 
ba \ 
LS \ BI (1 = L(eo44 —— 24) | {é of ee Cae > Ze 
Proof. We shall show by induction that 
k—-1 k—1 ) 
e+ LD altar —%) S| TT + heal le + Lesa} 
for k—1=7, r+1,...,¢; this will prove the lemma, since 1+L(z,,,—2,)=21 


10r 6 =1, 7=-41,..., t, and eqabe (era 0. 


For k—1=, the inequality follows from the hypotheses and the fact that 
1+L(2z,.,—4,) 21 because 


€ - Lg, (241 =f, = {1 =tg L(z,4 1 %,)} te Se Lg, (CA 2) 


Proceeding by induction, suppose the inequality holds for k<t—1. By the 
inductive hypotheses and the hypotheses of the theorem, 


k-1 
gies é = L 28 (Coy oe 2) 


But then : 
k k-1 
a L 2 Sx (2,5 2) = fe = L 2 Bae (Zat4 i 2,)} a Lg, (2p41 a %,) , 


and by use of d) this last term becomes 


eee ))} fe Fe oe lesa — 2,)} {4 car L (2,41 — 25)} 


M+ Leal etlelu-ah Ged 


5. Definition and Some Properties of p, and 4q,, 
If (x, y)ER, , let 


ae (x, y) if <4 <U44; 

Da; 9) = oe (7,9) if «=i; and 7<k, or if x=A; 
as) if w=t, and j>k, or if x=0; 

x 
(5.1) 

oe (x, y) if <y<taias 

Jn (%,¥) = rapt ave. ANG aks OF 1 oi) = 03 
a1 


“m (x, y) leet and Yoh). oral 7 =A, 
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“a 


We shall show that , and q,, which are defined throughout %, satisfy the 
hypotheses of ARZELA’s theorem (see Appendix 2). As a first step consider 
(x, y) ER, , and (%,v) €R;,,. Then 


By the linearity of the partial derivatives, the term on the right of the last 
inequality is the maximum of the four expressions 


otu OT Uy, 

(a) ae (t;, te41) = aq ie bet WI, 
ou, C Uy, 

(b) ‘ 0 “ (t;, ty) ee (t;, t,) ) 


or Un ( 


GU, 
Sere) ie 


+ |e 


t;, tya4) =y 


By (3.1) and (3.2), the terms in the last absolute value sign equal 


thy tet. 
thes face n) dn + S(t) —2(G)| SS |F(.n)|4nSM-|t..—-4|SMa,. 
tk 
Hence 
Ot Uy, OF Uy Of Uy * n 
|S (by taaa) — SE" A)] |S, tee) — SG tea) | + Mey, 


A similar result holds for (d). Thus we need only consider the terms (a) 
plus Mz,,, or (b) plus Ma,. Furthermore (a) arises from (b) by replacing k by 
k +4. Thus 


—" (x, y)|S max 


oe i (tis th) — EM (Gt) 
(5.2) On On Osea; —n|| ee *¢ eae wad 


+ M«u,, 


= imax |7( Xj, Vp) Tt(X;, V,)| + Ma, 
if xSy, x<y. The last equality does not hold if x>y or ¥>y or both, but the 
inequality between the left side and the right side of (5.2) holds for any (x,y) € Dia 
and (x, y) €9t;,. This can be shown as above by, for example, using (d) in place 
of (b) if Eee and Tie). 
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We now proceed to estimate the last expression in absolute value signs. 


| 20 (%5, Ya) — 2(%;, Ve)| Sl O71 —0 


Vk 
X; nan — J Flan) aa] 


ty 
S654. — ja] +| S/F) dn — fr Hee Olesaned 
by 


S2Mz, #28 (by A+ bal bleh E | Ope = Gj 4 


+] JPG ml a 


x7, ) —Flx;,m)| an], 


S 2M a, + | o741 — G1 ae i 


Gy) — flA(%;, 9), w(%;,0), »(%;,), (%;,), @(%7,)] + 
+ f[A(x;, Aeron: »1), %(%7,), @(%;,)] — 
— f{A(x;,), w(%;, 0), ¥(%;,7), @(%; 0), o(%;,n)]| an], 


te + op44— On| + MG —4] 4 


(45,7), M(%;, 0), V(X, ), (7,7), O(%;,)] — 


= HAC nh elm), 70) ela.) OC.) dn] + 
ae ES les — 9(%;,7)| an|| i ZS Lay.) — 2 (%;,1)| dn}. 


Let y, be chosen so that t;< yg<tg,,; then 


k—-1 


— 1(x;,%)| dn| = 2 (lor — te) |(%;, Vp) — 10(%;, Yp)]- 


Let y, be the terms in the brace on the right side of (5.3). Then, by use of Theorem 4 
with z replaced by ¢ and gy, by |2(x;, vs) —2(X;, ya)|, (5.3) can be written 


| x (x5, Vp) — 1(%;, Ve)| 


pas 
(5.4) Sy%,+ Ed pia — ts) | 2 (x5, Ve) — 7 (%;, Va) , 
= (Ud [1 + Lega — iI} {yy + L(ty41 — &) | 0 (%7, 97) — 2 (%;, y;)|t. 

By (3.2), (5.2) we have 


[Pn (%s y) — p,(%, y)| = ce (x, y) es OUn 


< max |2(;, y,) —(%;,¥,)| +Ma,. 


~0<k<n—-1 


(5.5) 
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We can now prove the following: 
Theorem 5. Given e>0, there exists 6>0, N2O such that if (x, y)ER,», 
(%, 7) CR; , and |x —x|<d, n2N, then |Pn(%, 9) —Pnl%, ¥)| Se. 


Proof. Choose e’>0. We will first show that y, and |2(x;, ;) — 1%, ¥;)| 
can be made less than e’ by taking » sufficiently large and |t;—4¢,| sufficiently 
small. 


| 1 (%7, V7) — 7%), V5) | [Op = Ga Fires ,4)| an =F |F (xn) an 


|p 0G ae praia + Mt, 


1— 4 
and clearly we can find 6’>0, N’20 such that the right side is <e’ whenever 
n= N’ and |t;—t,|<0’. 


We next consider y,. We have 


<|o,—o,|+2Ma,+M[a,+|5—4 


t% 


2 |o(%;.) — n)|an| S 


[seen ae— J FEM ree te 


> Tia , 


ty 


<M |t- || — 4] + | ta — Toal- 


Thus we see that we can find 6° >0, N’’ =O such that whenever n2N”, 


(a ll dn| S Be 
(b) aie 0;41| =e’, 
(c) M: |G —t4,| S50 


Consider now 
(e) [ALA (am), (a5, 09), » (45, 9), 20(%5,), 0 (45.0) — 
— FLA (% 1), eX), ¥ (%5), 7 (%5,), @(%5,0) I] 
We have already noted that z and @ are uniformly bounded for (x, y) ©. Also 
ylx — ‘4 = 2 
|y(%, 9)| S max, |un(,)| SM {4? + 24}. 


Since f(x, y, 2, ,q) is uniformly continuous in all five variables for (x, vy) ER 
and z, p,q bounded, we see that we can also make (e) Se’/4.4 for |x;—x,| <0" 
and n=N" by suitable choice of 6’>0 and N’>0. Thus y,Se’. Hence by 
(5.4) and (5.5), 


(Petry) Phan ys I [ee ibe {A+ Lys — )} + May, 
Se(1+La,exp{LA}+ Ma, 


Cauchy Problem for the Non-Linear Hyperbolic Equation 305 


whenever »2max {N’, N”}, |t;—t,|<min {6’, 6’. Since e’ can be chosen 
arbitrarily small, this proves our theorem. 

We are now ready to show that the sequence {#,! satisfies the hypotheses of 
ARZELA’S theorem. 


Theorem 6. The sequence {p,\ is uniformly bounded. 


Proof. |p, (x, y)|<f|F(x,y) |dy| + |2(4)| SM |x —y| +M<M(A +1), 
g.e.d: 2 


Theorem 7. Given e>0, there exists 6>0, N=O such that 


| Pn (%,) ae Py, (%, ¥)| <é 
whenever |x —x|<6, |y—y|<6, n=N. 
Proof. Let (x, y)ER,, (9) ENR; x. Then 


| Dn (%, 9) — >, (%,y)| = | Pn (%, 9) — Py, (XV) =e lpn (%, 9) == Pa (Oey): 


Now by Theorem 5, given ¢’>0, there exist 6’>0 and N’=0O such that if 
n=N’, |x—x|<6’, then |p, (x,y) —,(x,¥)|<e’. Hence 


Pn @H) — Pale Se +|f Fle 0) dy — J F(x) dn] 
Se+M\y—y|Se'+Mo, 


from which the conclusion of the theorem is immediate. 


We can now apply ARZELA’s theorem and conclude that there is a subsequence 
{by} of {p,} which converges uniformly to a continuous function # in the region 9. 

By symmetry it is clear that everything we have done for p, goes through 
for g,,, and hence there is a subsequence {g,-} of {q,,} which converges uniformly 
to a continuous function g in the region . 

If we first choose the subsequence {q,}, and then choose {f,,} to be a sub- 
sequence of {f,-}, an finally choose {u,,} to be a subsequence of {u,,} then we 
see that 


(5.7) lim 4,,= 4%, lim £,,=?, lim q,,= 4. 
6. Convergence to a Solution 
We wish to show that , the limit of the sequence {w,}, is a solution of the 
Cauchy problem. 
We have already shown that {w,}, {p,}, and {g,} are uniformly bounded for 
(x, y)ER, ve, there exists Z, P, Q, all 20 and independent of m, such that 
|un(% MSZ, |Pr(% WISP, |on(% y)[SQ. Also f(x,y, 2, 6,9) is uniformly 
continuous for (x, y) ER, |z| SZ, |p|SP, |g| SQ. For all n20 and for (&, 7) CR 
we have 
m(E,n) = Pa lA(E,n), u(E,n)], 
0 (6.9) =a ACE, m), u (Em) I, 
v(E,n) = u,[A(E, 7), w(E,m) | 


(the z, 0,v, 4 and w should also carry a subscript 2, but we omit it). 
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Also, by renumbering the elements of our sub-sequences {up (Dy 4dr te we 


can write 


(6.4) lim u,(x,y)=«(x,9), lim ,(%, 9) =p(%y¥), img, (x, ¥) =49(% 9), 


the convergence being uniform for (x, y) CM. 
Consider now 


ules) — {= fan fil8.n u En), POE mo 9(G-m)) 48 +fo6) aE+fx(8) a8} 
Sian Bete u(E,n), p(n, (&)] — 
—f{A(E,n), me (E19), Un (ACE, 9), WE, 0), Bu (ACE, 0) ol, 1), In (ACE, 1), oe(E0)) J} 
+ {2 — ods + f(T) — ce} ae. 


(6.2) 


Let us write A, u for A(E, 7), u(&, 7). Then 
£18, 9, #(E,), bE, 9), (E19) ] — FLA Ms Mn (As) Pn (As Ms In (A, 14) 
<|f[En, w(E,n), b(E 0), 9(E)] — tA wu), P(A), a(A Mw) + 
+ [FLA me 0 (A, ), BA, ), (A, H)] — FLA, Me Un (As He)» Pn (As tt) Gn (A, 12) J] 
Now by (6.1) and the uniform continuity of /, given e>0, there exists N’=0 
such that if ~= N’, then 
|7[A, ou, (A, oe), D(A, We), 9A, )] — FTA oe Mn (As) PAs) In (A) ]| Sse 
for all (A(E, ), w(&,)) ER 
Moreover from the definition of A, u it is seen that 


im A(E, 4) =§, lim w(E,9) =m, 


the convergence being uniform for (€,7)CR. Also w,, and qg are uniformly 
continuous in ft. Hence there exists an N’’>N’=0 such that if n=N” then 


[71é,9, #(E,7), B(&,n),9(€,9)] — fA, mw, w (Au), P(A, mM), 9(Aw)]| SHe. 


Again, since 2’ and T converge uniformly to o and Tt respectively, there exists 
N2N”2N'2Z0 such that if n2N,,|2(&) —o(€)|Se, | T(6) —1(6|Se for all p 
such that OS[€ <A. Hence for n=N we have by (6.2) 


U,(%, Y) — i Fan SHE, num), P(E, ), 9(E,) | dE + 


+ folgas + fre as\| <e{4* +24}, 
0 
for any (x,y) Ct. This shows that 


Gay vee — f dn fF18.n.u(8.n), pln). aem)ds +Jole) ae + frag 
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Hence we have 


_ (x,y) = filen u(x,7), P(x, n), 9 (x,y) ]dn + a(x), 


x 


2 (wey) = — ffl, y, wl.) bE»). g(E NEE Lely), 


oy : 
Cru Cu 


bx by Y= Fy aq % %) = flay, w(*, 9), B(x 9), 9(,9)]. 


All that remains to be done is to show that du/déx =, du/éy—=gq. To this 
end consider the following, where 2 =A(x,) and uw =(x,y): 


Pals 9) 4 f flan u(x, n), bm), ale ml dn +otx)h 


(6.4) = (FLA, os Mn (As MH), P(A 4)» On (As) ] — 


— f[%,n, u(x, ), P(x, 7), q(x, 7) |} dn + 2(x) — a(x). 


In the same way as for (6.2) we can show that given e>0, there exists an 
N’>0 such that whenever > WN’, the integrand on the right side of (6.4) is in 
absolute value less than e. Also if t,<*<t,,,, we have |2(x) —o(x)|=|o,,,—0(x)|, 
and since o is uniformly continuous, this will be less than ¢ for »=N’’ where 
N"’=0 is independent of x. Combining these results, we see that 


p(x, v) =f fle, n, w(x, 0), P(x, 7), 9(%,),] dn +(x) 


a) 
+ sia (x,y). 


x 


Qa 


Since g can be treated in an analogous manner, we have proved that a sub- 
sequence of “approximate solutions’, {w,,}, converges to a solution of CAUCHY’s 
problem. 


Section II. The Mixed Boundary Value Problem 


7. Statement of the Problem 
Section II deals with the mixed boundary value problem, the precise formula- 
tion of which is given in the following theorem. 


Theorem 8. Jf f satisfies conditions a), b), c), d) of Theorem 1, and if 
and wy are functions with together with their first derivatives are continuous real- 
valued functions of the real variable € for OSESA, and —y(0) =y(0), then there 
zs at least one u such that for OS x, ySA: 

(i) ws a function of x and y; 

(ii) « and its partial derivatives u,, Uy and U,,=U,, are continuous for OS x, 

V=A, 

(iii) 4, (%, ¥) =f[%, ¥, U(X, ¥), 4, (%, ¥), My(%, ¥) I]; 

(iv) «(0, 9) = p(y); 

(v) u(x, x) =y(x). 
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The solution u(x, y), which exists by Theorem 8, is not necessarily unique, 
as is shown by the following example. Let f(x, y, 4, £, g) =—4(7y?—4*y) | «|, 
Upy=—4{7y2—4xy}|ul*, u(x,0)=0, u(x, x)=0. Then u(x,y)=0 and 
u(x, y) =(v2—xy)4 are both solutions. However, as in the case of Theorem 1, 
if f(x, y, 2, p,q) satisfies a Lipschitz condition in z, in addition to the other 
hypotheses of Theorem 8, then the solution is known to be unique. This is proved 
in [5] for the linear case, but the proof is the same in the non-linear case. 


8. The Approximate Solution 
If / is a function of x, y alone, and is independent of z, #, g then a solution 
of the mixed problem is given by 


y % 


(8.1) w(x,y) =f dn f En) as + vm v(x) + ey). 


That this is so is verified as follows. Clearly u(0, y) =@(y), and u(x, x) =y(x). 
Also 


Hs(%9) =J Hn) dn — f Ea) as + ye) — 9), 


x 


(8.2) uy (x, ¥) = JS (Ey) dé + o'(y), 


0 
Ury(%,¥) = Uy(%Y) =H(%¥), qed. 

As in Section I, let ¢ be replaced by a function which is constant on each 
vj. Also let y and @ be replaced by functions which are continuous but with 
piecewise continuous first derivatives, 7.e., the first derivatives may have dis- 
continuities at E=t;,7=0,1,...,2—1. Then, again as in Section I, formula 
(8.1) would still give “solutions” to the mixed problem if (x, y)CR”,. 


We now proceed as before, introducing the “‘approximate solution” as 


Dh x 
(8.3) u, (x,y) = J dn f F(&,; x,y) dé +P, (x) — G, (x) + @,(y), 
x 0 
where F,, %,, and @, are defined by: 
Y; at p(t) ; Yj =y(i); 


(rent if 7 +k, R+0: 
Uy p= Q; if k==.0; 
i“ ih eae 
(i. if k>7 and (x, y) E|Rj i 
DOS) rare 
lujnir if RSJ and (x, y)E[R,, 
e if kj and (x, y)e[R,, 
oo if kj and (x, y)€[R,, 
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tkti—* 


P(x) Sm |; 2ae (Wpia aa Yr) if be <a 
Vn if ie 


tessa 


P(x) -| oy (aay (Peri— 9x) 1h  S¥ <tag; 
Pn if 16 Sh 


G4y—h 


iff<k and (x, 9) [Rj a5 
m(%,Y) = eS, 


Uti k+1_ “j,k41 rae ist 
eae ui 72k and (x, V) [Ry g- 


o(x, y)= Eire Af feco pecan (x, YER, a or if7=2k and (x, y) ER, y. 


teti—th 


O° if (2yER,, and (E,n)ER, 


Fin yi&m) =| ps 
f[A(x, y), A(x, y), v(x, y), BU vy), o(x, y)] otherwise. 


F*(x,y) =f[A(%, 9), (4, 9), ¥(% ¥), 2(%, 9), O(%, ¥)]. 


We next derive the analogue of equation (2.4). If7<k and (x, y)€ 
by (8.3) 


M,, , then 
y x 
u(x, y) =f dn f Faé + P(x\ — O(x) + Dy) 
x 0 


y x 
=fdnfPFdé + 
5 bj 


tk 


y tj 
+ fan [Pas +¥%—2,+ 0) + 


th x 


+ fdn Jl Fdé + P(x) —O(x) +O, 


t 


k bj 
a CE eB ae 


= fin (% t;) (y ty) u(t;, y) se u(x, tn) — U4 pr: 


We should get an analogous expression if (w, y)€ HR, and 7ZR. 


9. Some Properties of u,, 
If (x, y)ER,,, then by (8.3), 


x 


V, 
OUn (x, y)= it FG x) dy ip (ES x34, vide + P(x) — O (2), 


Ox ; 
(ioe , 
Stes (x, y= f FE, yi 4,9) dé + Oy). 
s 0 


Arch. Rational Mech. Anal., Vol. 3 25 


370 JAMES CONLAN: 


From the definition of F, ®, Y, we have for (x, y)CR,, 


otu = OUn CO Uy 
pak ed bene ie Naser ad Ua) 
(9.2) Ot Up OUy (x,y) = O Uy ie 238 
rn (x, y) = Sie (x, 9) = GO" (teas 
We also have 
a ee = 4, i (Pti41) bas” P(é;) a P(t, 41) a D(t;)} 27 
sal nee “iy 
te ttn 
sas i, =a Jan | ae pees 
‘Tere! 
ths 0 
(9.3) a m 
fy { fen fas fan fas} me 
bat 

Viti Vi Pj+1— Pj __ Un (x,, t,) 

tj 4—b; tay—tj ame nth 
where, as in Section I, t;< %;<4,41. 

Similarly we find that if («, y) CR, , 
Uj k+i—Uj,k _ OUn (4 
thigh ay (tj. 9%) 


where t,<. y<tjiy- 
As in Section I we choose a sequence of subdivisions of %t such that 
lim a, = lim max |t,—#4,|=0. 
n—> 00 n—>0o 0Sj<n—1! / oe 
Since ¢, gy’, y, y’ are continuous in the closed interval [0, 4], these functions 
are bounded in absolute value there. Let MW be a common bound of these four 


functions and of the absolute value of f. Then M will also be a common bound 
of |F |, |], |¥"|, |®|, |®’|. 


We now prove the analogues of Theorems 2 and 3 of Section I. 


Theorem 9. The elements of the sequence {u,\ are uniformly bounded in KR. 
Proof. By (8.3), u,(%, vy) SM{A+3}, q.e.d. 
Theorem 10. The sequence {u,} is an equicontinuous family in R. 


Proof. As in Theorem 3, assume that »<x, y<y. Further let us assume 
that *<y,x<¥y. 


| %n (#7) — tn (x, ¥)| S| PCH) — P(x)| + | B(x) — G(x) + | Py) — Py)| + 


ES 1 ae 
x 0 
S |P%) — P(x)| + |O(%) — O(x)| + | Py) — P(y)| + 
+ fa nf ag — fan fas FE, m| +20, M. 
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Consider the difference of the following double integrals: 


y x y x 
Jan fas sling) as| ONG n) 
min (%, y) 


{ fanpage’ | anfas) Pm 


(¥, y) x 


<MA|y—y|+MA|x —x|+MA|z%—x| 


<MA{\y—y| +2] — a1}. 


x 


‘|B E, m)| 4& 


Now given ¢>0, a exists 60> 0, ae such that if Pa best? ly—y|<6, 
n>WN, then |Y(x) — eee é, he (x)|<$e, 2a,M<te, MA {\y—y|+ 
2|%—a|}<te. ton | u,,(%, vy) —u, (x, aie queid: 


We see by (9.1) that 


au, 


aye y )| SM {|x —9| +|*|}+2M<2M{A +1}, 
and 
OUn 
|S (2,9) SM fA +4}. 


Hence z(x, y) and o(x, y) are uniformly bounded for (x, y) ER. 
The remarks following Theorem 3 in Section I hold unchanged for our present 
problem. 
10. Definition and Some Properties of p,, and q, 
If (x,y) ER, , let 


OUy, ; 
i eh ig et as 
OT Un . 
ie y) if x=t 
OUy . ; 
ay (x,y) if t.<Y <ty13 
/ Ot Uy, : ¢ ; 
(10.1’) In(%¥) = 45 (%y¥) iby=t and 7<k,; 
Se (x,y) iw y=h4, and 72k. 


As in Section I, we shall show that #,, and q,,, which are defined throughout %, 
satisfy the conditions of ARZELA’s theorem. Assume that (x, y)CR;, and 
(x,y)€R;,. We proceed exactly as in Section I, obtaining 


Ot Uy, OT Uy, 
| < (¢;, ty +1) a wenn (5, tr) 
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Using (9.1), (9.2) to estimate the last absolute value, we have 


(7 — f\Peman—{f— FE) dy sf iP n)| dn <Mo,. 


Continuing as in Section I we find 


OUn ( 


OUn (= \ 
(10.2) | in (%, 5) — 5 


an x, y)|s = "max | = (%;, Ve) — (Xz, Vx)| +Mz«,, 


O0sksn-1 


where, as usual, t;< %;<tj415 Gy< Vp <Aei- 
Proceeding in the same way, we could show that if (x, y)€ 2, , and (%, y) CN; - 
then 
, OUn ps = as OU, a 
(10.2") | St (2,59) = 2 (#,9)| s 
We next estimate the right side of (10.2’). 


| o(x;, Ve) 4 o(%;; ye) + Ma, 


ma 
Gotene ui 


|0(%;, vz) — o(%;, ¥a)| 


gins Vz) db — [FE Yi 4) Vr) agé 


a) 
<= |@'(y,) —O'(y,)| + ‘ 


S|'(yx) — O'(y,)| + 2M, ate |F*(E, ve) — F*(E, y4)| 26 


< {|©o%) — (| + 2Ma,+ 
(10.3’) 


SEG Va)» M(E, Ya), V(E, Ya), (EF, Ye), O(F, Ya) ] — 
— FACE, Ye), CE, Ya) ¥(E, Ya) 1(E, Ya), O(E, Ya) LAE + 
= LS |, yn) — @(E, ¥p)| as\ “ty L Slots yx) — o(8, ¥)| 46. 
As in Section I, we let y, be the term in braces, and we get 
|o(%,, Ye) — Q(x; Ya)| 


(10.4') aie bs {Ellas 4) [Olle Ye) — 0 (tes HL} 


< {TT 1+ Lle—4)]} 1+ Let) Lette. 99) — oll, wll} 


Theorem 11. Given e>0, there exists 0>0, N=O such that if (x, VER, p 
(%, vy) ER, 5 and 
lIy—y|<6, 2N, then |¢,(%,¥) —¢,(%.9)|S 
Proof. We have 


| On (X, ¥) Sadi x, y) Js < max |o(x;, VR) — 0(%;, Ve) | EG, 


0SjSn-1 


={ iT [4 + L(t py 7 A} fntLa— to) | 0 (fo, Va) — e(%, Vx) |}- 
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As in the proof of Theorem 5, we show that given e’>0, we can make y, 
and |o(t, Vz) —@ (to, ¥;)| less than «’ by taking sufficiently large and |¢; —¢,| 
sufficiently small. In fact 


| oto, Ys) — (to, Ya)| =|P'(va) — O(y,)| 


Prt+i— PR Pkrt+i— Pk 


thii—te te+1—"e 


and by the mean value theorem this is 


5% tts +8) — 5° (4 +8)! 


where 0<#< (tg41—te), O<O<(th14—4&). But q’ is uniformly continuous in 
[0, A]; hence we can take m large enough so that 


OD jo i dq e’ 
Ee (tr ) dy (tg) Se 
dy dp e’ 


Hence 


Jello. 98) — ello, wl S| 4% (i) — 45 


and we can take |¢; —¢,| small enough so that the right side is Se’. 
Next consider y,. We treat |®’(yz) —@’(y,)| as above. Also, 
Vk 


Vk 
SF(E, 7; &, yg) dn — FFE: 8 9) dn| df 


& 


Slate, yi) =m(E, | des f 


tj 
<2Ma,|t;| +f 
0 


Vk 

[\FAE | dn] ae 
Vk 
<2Ma,A+M|y,_—y,| -|¢,| 
<MA{2a, +|vs—%l}- 


We now proceed as in the proof of Theorem 5, getting y,<e’ for suitable choice 
of m and |t; —¢,|. 


Hence 


| an.(%,.9)) — Ga (*, ¥)| = max | o(%;, Vs) — 0 (%;, ¥s)| + Ma, 


~ 0sjsn—-1 


Ma, +{ TM + ite —t)lh eft + Ls — 6) 


<Ma,+e«'{1+Lza,}exp{L A} 


for n sufficiently large and | t; — t,| sufficiently small, and this proves the theorem. 
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We next estimate the right side of (10.2). 


| 0 (%;, Ve) — (45, Vp) 
<S [W"(x;) — P'"(x5)| + 19'(x)) — '9)| 


oe eee 


BAG riko 
0 
eae A 4%;)| + 


a nd + FP, pat 


ai |E*(x;,m| 4m — \Fl(.n)lan + 


(x3, + FAC Nan] + 


ae | fF.) F |F¥(x;,m)]} dn] +| UP am Si ro’ 


hj 
+ |S Fee; t) 4 
0 


(Et, yas] +], 


(&, &) as| 
SP" a) — P(x) + [P') — Op] + 
(10.3) + M{|tui—4| +|441—-4| + 2/4 — foal + 
2|Vp t,| 2| tray — Yel + |4;—4| + |= Gl} + 


ty 
+ | SPH a5.) dy 
by 


(x;,m) dn] + 


yy 
te | LER E, Glee 
0 


(35.10) — F(x.) dn 


<{|P") — "'(x,)| + |G'(x;) — O'(x,)| + M[10e, + 2|4—4]] 4 


+ S|FME.4) PME ylagh + 
0 


LIP*ep — F*(x;,n)| an) 
< '(x,) —W"(x,)| + |®' (x) — O'(x)| + M [100 + 2|&—4|] + 


4 
Ta [F*(E gy —FA(E, t)| dé 


+] SUA Gp) selon), 5.2), 201%.) 0-7] — 


— f{Alxj.n), ele). 0), ¥( 1), 22%), 0 (35, I) a) + 


+ [EFetsnn — ets an]h + [Ef xm) — xm an] 
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As in Section I, we let y, be the terms in the last brace, and we get 
| 2 (x;, 0 ie Te (%;, Vn) 


h-1 
(10.4) Se E ds ors re ts) | (>, Vp) To 76 (%; , Va) 


= {i lp ae AC) ime tp))} {yr + L(t 14 — &) | a (%;, V;) — 10(x;, yz). 
=, 
Theorem 12. Given e>0, there exist O0>0 and N>0 such that if (x, y)ER, », 
(x, V)ER,, and |x—x|<6, n>N, then |p,(%, y) —, (x, y)|Se. 
Pioof. 


Sees eat Vr) am me (%;, Vp) ae Mu, 


= (L uae Lite i4 =F in)}} {ny ae L(t+4 = t;) | x(x;, 7) = 1 (%;, y;)|} ee 
Sl 
We proceed as with Theorems 5 and 11. Choose ¢’>0. 


| 2 (Xx;, Yj) ie Te (%;, 4;)| = 


NF XG 
<2Ma,+ J |F(é, y;)dé|S2Ma, + M|x,— x; | 
x 


and this is <e’ for » sufficiently large and |x;— x,| sufficiently small. 
Now there exists 6’>0, N’>0 such that whenever n> WN’, 


[Y"(34) —'(%,)|< zee’, | O'(%;) — O'(x))| <tee’, 


x;—x,;|< 6’, then 


which is shown in the same way as in the proof of Theorem 11; and 


te 
E| {lols 0) — e(%.m)| dn) = ree’, 
“if 


which is shown as in the proof of Theorem 5; and the expression (e) in the proof 


of Theorem 5 is S Sy e’, and similarly if the 7 and x; are interchanged; and 


M[10%, + 2|%—4]] < zee’. 
Hence 


ij 
neste +f \FG 4) — Pe dts ge + te’ 
0 


vy) bj 
+L J |x, 4) — 26,4) dE +L J lo(€,4) — 06, 4)| 46. 
0 0 
As in the proof of Theorem 11, the first integral on the right is 
< LMA (2a, +|¥s— Yal}- 


By Theorem 11, the integrand in the last expression can be made < ior e’ for 


proper choice of |¢;—#,|. 
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Combining all this, we see that y,<e’ whenever n>N’, |;—#,|<06’. Hence 


|Pn(%,¥) — bal(% ¥)|S max | 0(%;, Vg) — 1 (%;, Y—)| + Ma, 


~ 0SkSn-1 
< {11 + Lp —o)] e+ LG — 9) + Mas 
si Ra + {1+ La,}exp{L A} 
or n sufficiently large and |t;—¢,| sufficiently small, which proves the theorem. 


Theorem 13. The sequences {p,\ and {q,} are uniformly bounded. 


f(x, | an] + SLPME, 2) dE + 12") + 0) 


Proof. |, (x, 9)| < 


<M {|x —y|+ |x| +2}<2M {A +14}; 
and 


lan(4, ”)| SSIF*E,y)| dé + |O'(y)| SM{A+1}, ged. 
0 
Theorem 14. Given e>0, there exist 6>0O and N>O such that 
[Pn(¥%¥) —Pul%M1<e, |an(®V) —Anl¥% | <E 
whenever |%¥ —x| <4, |y—y|<0, n>N. 
Proof. Let (x, y) ER, ,, (x,y) CR g- Then 
|Pn(%¥) — Pal, ¥)| S| Pn (%,¥) — Pn (%,¥)| + | Pn (¥ 9) — Pa(%, ¥)|- 
By Theorem 12, given e’>0, there exist 6’>0 and N’>0 such that ifn>N’, 
|x —x|<06’ then |,,(%, vy) —p, (x, y)|<e’. Hence 
y y 
|Pn(¥,9) — Pa(¥,9)| Se’ +| fF*(x,n) dn — f F*(x,n) dyn) + 2Ma,, 
x x 
Se’+M|y —y| +2Ma,, 
which proves the theorem for {#,,}. 
Similarly 
| Qn (%.9) — In(%,¥)| S| On (¥, 9) — On (%,¥)| + | Qn (2%, ¥) — Gn (%, 9) 


We apply Theorem 11 to show that the first term on the right is <e’ for 
n> N’', |y—y|<6'. The second term is 


S| J F*(é,y) d&— f F*(E, y) dé] +2Ma, <M {|x — x| +20,}, 
0 0 


which essentially completes the proof. 


As in Section I, ArzELA’s theorem now guarantees the existence of sub- 
sequences such that 


(10.5) iim 4,,=4, lim p,=p, lim ¢,,=4, 


the convergence being uniform in %t, and the limit functions being uniformly 
continuous in i. 
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11. Convergence to a Solution 
We proceed as in § 6. However, instead of (6.2) we have 


alt —{ fan film w(E.n), bE n), a6 148 +y2) — 9%) + 90) 


(11.1) = — fan f a8 (/18,n, 0&2), PLE) aM —F,(E,n; x, y)}+ 
+ (x) — p(x) — ©, (*) + o(%) + ©, (y) — oly). 


With the same notation as in Section I, 


[flé,n, #(&,n), me n),9(€,)] —F,(é,n; *, 9)| 
S|f(é., u(€,7), 0(6,7), (6, 0)] —F LA Me Un (As), Pn(A, H)s In (Asm) | + Moe, 


and the term in absolute value signs on the right can be treated as in Section I, 
giving the result that for any «>0, there exists N’>0 such that if »=N’, then 


y a 
Jan J a&it[E 1, (En), P(E 1), 96,1) — FE m: x9) SeA?+ May. 
Moreover, if t, 5 *<t,4,, then 


|'¥.(%) — p(x)| = 


teti— x 
Le41—tr 


? 


Prt (74 =| (Vara — Ya) — P(A) 


and by the uniform continuity of y, there exists an N’’>0 such that if n=N”, 
then | p(x) —Yriil<Be, |Yari—Yal<de- Hence for n>N” |¥, (x) —y(x)|<e, 
and similarly for |®, («) —@(x)|. 

Therefore for n= N =max{N’, N”} 


U,, (X,Y) —{ fan free n, u(&,n), b(&, 7), 9(&,)] dE + p(x) — v(*) +o(y)h| 


< e{A?+3}+ Ma, 
for any (x, v)€ 8, which shows that 


(41.2) (x,y) =fanf tlm u(é,), P(&, 7), 9(&,n)] dé + p(x) — p(x) + v(y). 


Therefore 


oe (x, y)= f fle n, w(x, 0), b(% 0), ¢(% 0) ] an — 


— fil, x, w(E. x), DUE), 916.) 48+ ya) — 90); 


ay =f HB 9) PE 9). ENE + 9) 
fe 2H = I(x, 9, u(y), POY) 19): 


We still must show that du/dx =, du/dy =q, but the proof of this is essen- 
tially the same as in Section I. 
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Appendix 1. Numerical applications 
We apply the “approximate solutions” given by formulae (2.3) and (8.3) to 
the equation of LIouvILLE [4]: 
(a) Uny =exp{2Au}, A = constant. 


The general solution of this equation is 
ee In| 8 (vy) h'(#) } 
2h [Ag(y) +h (*)]? 


where g and / are arbitrary functions having suitable differentiability properties, 
and such that u, u,, and u, exist. This can be readily checked since 


(b) “= 


_ 1 fax) n(x) 
(c) Me 25a) 7? tetas? 
ee te fer) Rial a, \ 
¢) "21 tO) ely) +h(a)S” 
Wry eg gE) =exp{2Au}. 


"* They) +h (a)? 


The numerical process used in the proof of Theorem 1 will now be applied 
to compute numerically the particular solution of equation (a) which is obtained 
by taking A=1, h(x) =exp x and g(y)=sin y. The fact that the solution is 
actually known in this case will serve as a check on the computations. The values 
of o(&) =u, (&, é) and of 7(é) =u, (€, &) can be computed from (c) and (d) respec- 
tively. For our example we let R be the square with vertices at (0, 0), (0, 2/4), 
(7/4, 0) and (z/4, 2/4), and we subdivide ® into nu? squares. 


By (2.4) 
1 
(e) 4 5= 2D (2G) + TE) Ga —4) 
i=0 
and for 7<k 
(f) Ups = (Gy = baa) bea = Fe) Tee NG po pe ee 


Using (e) and (f) to compute w(0, +2) we get: 


LOR Hie, u(0, 4%) = — 0.7097; 
for n=5, u(0,+2%) = — 0.7092; 
for m=10, (0,42) = — 0.7088; 
for m=20, (0,42) = — 0.7085 


On the other hand, from (b) we have to 4 decimal places (0, +) = — 0.7082. 

Next we apply the numerical process used in the proof of Theorem 8 to 
compute numerically the particular solution of equation (a) which is obtained 
by taking A=1, h(x) =exp x and g(y) =sin y-+2y. We can then obtain x (0, y) 
and u(x, x) from equation (b). We let be the square with vertices at (0, 0), 
ee (2, 0) and (27, 2%), and we subdivide Rt into n? squares. By (8.4), 
if 7< 


(g) Cb May joa) (t44 ae t;) (te41 — &) hier ae Up Uj +41,n — 4; R- 
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Using (g) to compute u(x, 27) we get: 


for n=2, u(m, 27) = — 1.533; 


for  =10, 


( ) 

for n=5, u(t, 2%) = — 1.473; 
w( ) 

for 7—= 20, u( ) 


On the other hand, from (b) we have to 4 significant figures u (2, 2) =—1.454. 


Appendix 2. Arzela’s theorem 


Arzela’s theorem. Let {f,} be a sequence of real-valued functions, defined on 
a compact subset & of a separable metric space, such that: 


(a) The sequence {f,(x)} 1s uniformly bounded in absolute value on 8. 


(b) Given any e>0, there exists N>0O, d>0 such that whenever n>N, 
|%1— %_| <0, %ER, x2ER then | fF, (%1) —f,(%2)|<e. Then there exists a subsequence 
{t,,} of {f,} which converges uniformly to a continuous function g defined on &. 


Proof. Since St is separable, there exists a countable set BC such that the 
closiire of 8 is GQ. Let Wie, ¥,..00}. Choose a subsequence of Un} which 
converges on B as follows. Assume we have chosen a subsequence {f,, ;} } } of {f,} 
which converges at %,, %,,...,%;. Consider the point x;,,. Since |/, ;(;+,)| is 
bounded there will be a subsequence {f, ;,,} of {/,,;} which converges at %, 
Xg, +++, %;, %4,- Consider a new sequence {f,, ,,} which is formed by a diagonal 
process (ron the subsequences {f,, 1}, {f,,2},-... This diagonal sequence con- 
verges on B, 1.¢., 

im fa, n(% ) = g(x;) 
for all x,€ B. 

We now show that {/, ,} converges uniformly on st. By hypothesis (b) for 
any ¢>0, there exists 6>0 and N>0 such that if x, yER, |x—y|<dandn>N 
then 


(1) tna) ag to. Ps yyl<4 é. 
Also, since & is compact, for any 6>0 there exists a finite set {y,,..., y,}€B 
such that for any xG€R we have |x —y;|<0 for some 7 such that 157<p. We 
can take N large enough so that in addition to (1) we have for n, m>N 
Lie NO] arm (Ol)| owe 7=1,2,...,p- 
Then for some 7 (17) we have 
pase n 1, m(% x) | = Tegal, n(¥;)| Shee OA acs Pen OF Tallent AS aes) 
<= &, 


which proves the uniform convergence of {f, ,} to a function g defined on &. 
The continuity of g follows directly from hypothesis (b). 

AscoLi’s theorem is obtained from ARzELA’s theorem by specializing hypo- 
thesis (b) to read: ‘“The functions {f,} are equicontinuous on St.” 


380 James ConLtan: Cauchy Problem for the Non-Linear Hyperbolic Equation 


References 


[1] Leenery, P.: On the existence of not necessarily unique solutions of a classical 
hyperbolic boundary value problem for non-linear second order partial dif- 
ferential equations in two independent variables. Ph. D. thesis, Brown Uni- 
versity, June 1950. 

[2] Harrman, P., & A. WintNER: On hyperbolic partial differential equations. Amer. 

J. Math. 74, 834—864 (1952). 

[3] Diaz, J. B.: Onan analogue of the Euler-Cauchy polygon method for the numerical 

solution of Uyy=F (x, YU, Uz, Uy). Arch. Rational Mech. Anal. 1, No. 4, 

357 —390 (1958). 

[4] Jorvan, C.: Cours d’Analyse, 3rd edit., vol. 3, pp. 369—371. 1945. 

[5] Picarp, E.: Lecons sur quelques types simples d’équations aux dérivées partielles 
avec des applications a la physique mathématique. 


University of Maryland 
and 
U.S. Naval Ordnance Laboratory 


(Received February 19, 1959) 


at 


BVDITORTAL BOARD * 


S. BERGMAN 
Applied Mathematics and Statistics Laboratory 
Stanford University, California 


R. BERKER 


Technical University 
Istanbul 


L. CESARI 
Research Institute for Advanced Study 
Baltimore, Maryland 
4 Le COLLALZ 
Institut fiir Angewandte Mathematik 
Universitat Hamburg 
A. ERDELYI 


California Institute of Technology 
Pasadena, California 


J. L. ERICKSEN 
The Johns Hopkins University 
Baltimore, Maryland 
G. FICHERA 


Istituto Matematico 
Universita di Roma 


R. FINN 


California Institute of Technology 
Pasadena, California 


IRMGARD FLUGGE-LOTZ 


Stanford University 
California 


HILDA. GEIRINGER 


Harvard University 
Cambridge, Massachusetts 


i H. GORTLER 


pi Rerent fiir Angewandte Mathematik 
Universitat Freiburg i. Br. 


D. GRAFFI 


Istituto Matematico ,,Salvatore Pincherle“‘ 
Universita di Bologna 


A. E. GREEN 
King’s College 
Newcastle-upon-Tyne 
J. HADAMARD 
Institut de France 
Paris 


L. HORMANDER 


Department of Mathematics 
University of Stockholm 


M. KAC 


Cornell University 
Ithaca, New York 


H. LEWY 


University of California 
Berkeley, California 


A. LICHNEROWICZ 
Collége de France 


Paris 
GG; LIN 


Massachusetts Institute of Technology 
Cambridge, Massachusetts 


C. LOEWNER 


Applied Mathematics and Statistics Laboratory 
Stanford University, California 


W. MAGNUS 


Institute of Mathematical Sciences 
New York University, New York City 


G. C. McVITTIE 


University of Illinois Observatory 
Urbana, Illinois 


J. MEIXNER 


Institut fiir Theoretische Physik 
Technische Hochschule Aachen 


C. MIRANDA 


Istituto di Matematica 
Universita di Napoli 


C. B. MORREY 


University of California 
Berkeley, California 


C. MULLER 


Mathematisches Institut 
Technische Hochschule Aachen 


W. NOLL 


Carnegie Institute of Technology 
Pittsburgh, Pennsylvania 


A. OSTROWSKI 


Certenago-Montagnola / Lugano 


R. S. RIVLIN 
Division of Applied Mathematics 
Brown University, Providence, Rhode Island 
M. M. SCHIFFER 


Stanford University 
California 


J. SERRIN 


Institute of Technology 
University of Minnesota, Minneapolis, Minnesota 


E. STERNBERG 


Division of Applied Mathematics 
Brown University, Providence, Rhode Island 


R. TIMMAN 


Instituut voor Toegepaste Wiskunde 
Technische Hogeschool, Delft 


W. TOLLMIEN 
Max-Planck-Institut fiir Str6mungsforschung 
Gottingen 
A. TONOLO 
Seminario Matematico 
Universita di Padova 
R. A. TOUPIN 
Naval Research Laboratory 
Washington 25, D.C. 
CaTRUESDELE 
801 North College Avenue 
Bloomington, Indiana 
Ee VLE eA 


47, bd. A. Blanqui 
Paris XIII 


TLED IN SYACKS 


CONTENTS 


CoLeMaN, B.D., & W. Nott, On Certain Steady Flows of General Fluids 

Riviin, R. S., The Constitutive Equations for Certain Classes of De- 
formations SHE. Ws Ns} acon Sas a eel 

Tavs, A. H., On Circulation in Relativistic Hydrodynamics . 

Ostrowski, A. M., On the Convergence of the Rayleigh Quotient 


Iteration for the Computation of the Characteristic Roots and 


Vectors. III. (Generalized Rayleigh Quotient and Characteristic 
Roots with Linear Elementary Divisors) 

Ostrowski, A. M., On the Convergence of the Rayleigh ‘Giatient 
Iteration for the Computation of the Characteristic Roots and 
Vectors. IV (Generalized Rayleigh Quotient for Nonlinear Ele- 
mentary Divisors) Peni oe 

Weston, J. D., Characterizations of ees Transforms and Perfect 
Operators . 

ConLaN, J., The oie Problem and the Mixed eseinaaly Value 
Problem for a Non-Linear Hyperbolic Partial Differential Be 
in Two Independent Variables . 


Druck der Universitiitsdruckerei H. Sttirtz AG., Wiirzburg 
Printed in Germany 


289 


304 


342 


325 


» 341 


348 


1355 


